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1.  INTRODUCTION 


Gradient  models  are  analytic  representations  of  the  spatial  distributions  of  density, 
pressure,  temperature,  and  velocity  of  expanding  propellant  gases  between  the  gun  breech 
and  the  projectile  base.  In  a  thermodynamic  interior  ballistic  model,  only  mass  and  energy 
conservation  are  explicitly  required,  and,  therefore,  a  gradient  mode!  must  be  provided  to 
account  for  the  spatial  variation  in  gas  momentum.  Closure  is  achieved  by  including  a 
constitutive  law  for  the  gas  and  requiring  that,  at  any  point  in  time  during  the  interior  ballistic 
process,  an  integral  average  of  the  local  thermodynamic  properties  of  the  gas  equal  the 
corresponding  average  gas  properties  derived  from  the  constitutive  law  and  energy 
conservation  for  the  system  as  a  whole. 

In  general,  use  of  a  gradient  model  results  in  fixed  gradients  of  system  properties 
throughout  the  interior  ballistic  process  and,  thus,  a  fixed  ratio  of  projectile  base  pressure  to 
breech  pressure  (Corner  1950;  Vinti  and  Kravitz  1949).  Exceptions  arise  if  boundary 
conditions  and  ratio  of  propellant  charge  mass  to  projectile  mass  (i.e.,  charge  to  mass  ratio 
[C/M]),  vary  during  the  interior  ballistic  process,  as  in  the  Modified  Lagrange  gradient  model  of 
Morrison  and  Coffee  (1990),  or  if  multiphase  effects  are  included,  as  in  the  RGA  gradient 
model  of  Robbins,  Anderson,  and  Gough  (1990).  The  fixed  (or  constant)  gradient  character  of 
basic  gradient  models  is  a  result  of  the  exclusion  or  wave  phenomena,  such  that  changes  in 
the  physical  state  of  the  system  are  instantaneously  communicated  throughout  the  gas  in 
accordance  with  the  distribution  functions  of  gas  properties  for  a  particular  gradient  model. 
Thus,  gradient  models  entail  an  implicit  assumption  of  infinite  gas  sound  speed. 

The  development  of  gradient  models  has  been  closely  tied  to  the  ballistics  of  conventional 
guns,  and  such  models  do  not  appear  to  have  been  widely  applied  to  novel  propulsion 
concepts  such  as  the  light  gas  gun  (Seigel  1979).  However,  in  the  derivation  of  traditional 
gradient  models,  i.e.,  the  Lagrange  (Corner  1950)  and  the  Pidduck-Kent  (Vinti  and  Kravitz 
1949)  models,  there  are  no  assumptions  which  limit  their  application  to  conventional  guns  or 
to  gas  properties  characteristic  of  conventional  gun  propellants.  Therefore,  traditional  gradient 
models  should  be  applicable  in  thermodynamic  interior  ballistic  models  of  gas  dynamic 
propulsion  processes  irrespective  of  molecular  weight  or  sound  speed  of  the  driving  gas, 
subject  to  the  usual  consideration  of  expansion  ratio  and,  in  the  case  of  the  Lagrange 
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1950;  Vinti  and  Kravitz  1949).  Indeed,  the  accuracy  of  a  gradient 
configuration  should  improve  with  increasing  sound  speed  since  wave 
olumn  is  reduced.  As  a  result,  the  relaxation  of  wave  processes  which 
and  the  subsequent  establishment  of  a  nearly  constant  base-to¬ 
ld  occur  more  quickly,  thus,  at  a  lower  expansion  ratio  for  a  low 
[ising  the  range  of  validity  and  the  accuracy  of  the  gradient  model, 
pproximation  and  the  accuracy  of  ballistic  solutions  obtained 
Is  should  actually  improve  as  molecular  weight  is  decreased 


It  has  been  suggestvMthat  the  electrothermal  gun  concept  may  offer  a  practical  means  to 
employ  lew  molecular  weiljt  gases  (Juhasz  et  al.  1988).  In  the  case  of  an  inert  working  fluid, 
it  has  been  experimentally  Remonstrated  that  the  rate  of  conversion  of  liquid  to  gas,  as  well  as 
the  pressure,  density,  and  £  Serature  of  the  gas,  are  dependent  on  the  temporal  and  spatial 
distributions  of  the  plasma  s3|ce.  While  experimental  data  does  suggest  some  success  in 
achieving  ballistic  control  and  r|ueatabi!ity  with  inert  working  fluids,  the  data  does  not 
constitute  a  clear  demonstration  Ifc  ballistic  control.  In  the  case  of  energetic  working  fluids,  or 
propellants,  the  evidence  of  a  ci;  \correspondence  between  the  rate  of  conversion  of 
propellant  to  gas  and  the  characi^otics  of  the  plasma  source  is  less  conclusive,  particularly 
for  cases  in  which  the  ratio  of  propellant  chemical  energy  to  plasma  energy  is  greater  than 
about  2  to  4.  Although  experimental  data  is  somewhat  limited,  it  would  appear  that 
reasonable  ballistic  repeatability  may  be  achievable  for  ratios  of  propellant  chemical  energy  to 
plasma  energy  at  least  on  the  order  of  one. 


If,  however,  it  is  assumed  that  for  some  combinations  of  plasma  source  and  working  fluid 
(or  propellant),  gas  generation  can  be  precisely  controlled  by  the  plasma  source  such  that 
variations  in  the  rate  of  plasma  deposition  will  result  in  corresponding  variations  in  chemical 
energy  release,  then  repeatable  control  of  the  ballistic  process  may  also  be  assumed.  If  it  is 
further  assumed  that  such  permissible  combinations  of  plasma  source  and  working  fluid  are 
possible  for  a  sufficiently  wide  range  of  working  fluid  characteristics,  it  may  then  be  possible  to 
identify  unconventional  propellants,  which  upon  reaction  generate  low  (average)  molecular 
weight  gaseous  products.  In  comparison  with  the  combustion  products  of  conventional  solid 
propellants,  a  low  molecular  weight  gas  contains  a  greater  number  of  molecules  per  unit 
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mass,  the  ratio  of  the  number  of  moles  per  unit  mass  being  inversely  proportional  to  the  ratio 
of  molecular  weights  of  the  two  gases.  As  a  result,  for  a  given  gas  temperature,  the  low 
molecular  weight  gas  will  have  higher  specific  internal  energy  and  impetus,  as  well  as  higher 
speed  of  sound,  such  that  a  lower  mass  of  propelling  gas  is  required  to  achieve  a  given  initial 
gas  internal  energy.  Therefore,  for  fixed  initial  gas  internal  energy,  the  kinetic  energy  of  a  low 
molecular  weight  propelling  gas,  at  any  given  projectile  velocity,  represents  a  smaller  fraction 
of  the  total  chemical  energy  of  the  system,  resulting  in  a  reduced  projectile  base  to  breech 
pressure  ratio  {after  relaxation  of  wave  processes)  and  increased  thermodynamic  efficiency  of 
the  expansion  process  (i.e.,  higher  ballistic  efficiency  than  is  achievable  with  conventional 
propellants). 

In  recent  discussions  of  the  electrothermal-chemical  gun  concept,  it  has  been  suggested 
that  the  reduced  pressure  gradient  in  a  gun  using  a  low  molecular  weight  gas  is  a  result  of  the 
higher  sound  speed  of  the  gas,  which  permits  more  rapid  transit  of  pressure  waves  between 
the  breech  and  the  projectile  base.  Furthermore,  it  has  been  suggested  that  traditional 
gradient  models,  which  exclude  wave  phenomena,  should  not,  therefore,  be  used  in  interior 
ballistic  simulations  of  the  electrothermal  gun.  As  noted  earlier,  the  derivations  of  such 
gradient  models  do  not  introduce  limitations  on  the  properties  of  propelling  gases.  Thus,  it 
would  appear  that  their  application  should  not  be  limited  to  the  conventional  gun  case. 

The  objective  of  this  investigation  is  to  computationally  explore  the  range  of  applicability  of 
the  standard  gradient  models,  Lagrange  and  Pidduck  Kent,  in  the  simulation  of  interior  ballistic 
processes  involving  gases  of  varying  molecular  weight.  To  simplify  the  analysis,  the  Lagrange 
ballistic  problem  is  used  as  the  framework  for  simulations.  The  lumped  parameter  IBHVG2 
computer  model  (Anderson  and  Fickie  1987)  was  used  to  perform  interior  ballistic  simulations 
with  the  Lagrange  and  Pidduck-Kent  gradient  models.  The  one-dimensional  XNOVAKTC 
computer  model  (Gough  1990)  was  used  to  provide  wave  dynamic  solutions  as  a  baseline  for 
comparisons.  The  results  indicate  that  traditional  gradient  models  are  indeed  applicable  in 
interior  ballistic  simulations  of  gun  concepts  utilizing  low  molecular  weight  gases. 

2.  LAGRANGE  BALLISTIC  PROBLEM 

The  Lagrange  problem  is  an  idealization  of  the  interior  ballistic  process  in  a  gun.  It  is 
assumed  that  the  gun  is  a  right  circular  cylinder  (i.e.,  chamber  and  bore  diameter  are  equal) 
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closed  at  one  end  by  the  breech.  It  is  also  assumed  that  prior  to  the  start  of  projectile  motion, 
the  propellant  is  instantaneously  burned,  such  that  the  chamber  is  initially  filled  with  a  gas  of 
uniform  pressure,  density,  and  temperature.  The  gas  is  assumed  to  be  inviscid  and  heat  loss 
to  the  walls  is  neglected. 

Comer  (1950)  has  provided  a  historical  summary  of  efforts  to  solve  the  Lagrange  problem, 
including  the  development  of  the  Lagrange  and  Pidduck-Kent  gradient  models.  These  models 
are  addressed  in  greater  detail  in  subsequent  sections  and  in  Appendices  A  and  B.  The  first 
wave  dynamic  solution  of  the  Lagrange  problem  was  presented  by  Love  and  Pidduck  (1922). 
They  analyzed  the  Lagrange  problem  using  the  method  of  characteristics,  tracing  the 
rarefaction  wave  generated  by  acceleration  of  the  projectile,  as  it  moves  through  the  gas 
column  between  the  projectile  base  and  the  breech.  More  recently,  the  Lagrange  problem 
and  the  solution  of  Love  and  Pidduck  (1922)  have  been  used  as  a  test  case  for 
multidimensional  interior  ballistic  models  (Schmitt  and  Mann  1981;  Edelman,  private 
communication). 

The  pressure  distributions  behind  the  projectile  at  various  times  during  the  ballistic 
process,  as  obtained  by  Love  and  Pidduck  (1922),  are  shown  in  Figure  1.  The  initial  chamber 
pressure  is  6,333.3  kg/cm2  or  621 .09  MPa.  As  the  projectile  begins  to  accelerate,  a 
rarefaction  wave  is  generated  at  the  projectile  base  and  begins  to  move  toward  the  breech. 
The  pressure  gradient  between  the  projectile  base  and  the  wave  front  prior  to  the  first  wave 
reflection  from  the  breech,  and,  thus,  the  gas  kinetic  energy  distribution,  are  characteristic  of  a 
simple  wave.  The  wave  reflects  from  the  breech  as  a  rarefaction,  producing  a  nearly  flat 
pressure  distribution  and  a  base-to-breech  pressure  ratio  near  1.0  when  the  wave  reaches  the 
projectile  base.  This  reduced  pressure  gradient  results  in  reduced  gas  acceleration  and,  thus, 
a  reduced  rate  of  increase  of  gas  kinetic  energy.  Upon  reflection  of  the  rarefaction  wave  from 
the  projectile  base,  the  process  repeats;  however,  the  pressure  gradient  is  substantially 
reduced.  By  the  second  reflection  from  the  projectile  base,  the  wave  disturbance  has  relaxed 
significantly,  and,  as  can  be  seen  in  Figure  2,  the  ratio  of  base-to-breech  pressure  is 
approaching  a  constant  value  of  approximately  0.9. 
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Figure  1.  Pressure  Distribution  Behind  the  Projectile  at  Various  Stapes  of  Its  Motion,  in  the 
Example  Studied  by  Love  and  Pidduck  (1922). 
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Figure  2.  Time  Variation  of  Results,  In  the  Example  Studied  by  Love  and  Pidduck  (1922). 
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3.  LAGRANGE  GRADIENT  MODEL 


The  Lagrange  gradient  model  is  the  instantaneous  solution  of  the  continuity  and 
momentum  equations  for  a  gaseous  propelling  charge  in  the  Lagrange  gun,  with  the 
assumption  that  the  gas  density  is  uniform  between  the  breech  and  the  projectile  base  {see 
Appendix  A).  This  leads  to  a  gas  velocity  distribution  which  varies  linearly  with  position  along 
the  tube,  from  zero  at  the  breech  to  the  projectile  velocity  at  the  base  of  the  projectile.  The 
resulting  pressure  distribution  is  a  quadratic  function  of  distance  from  the  breech.  As  noted  by 
Comer  (1950),  the  term  'Lagrange  approximation'  is  applied  to  both  the  assumption  of  zero 
gas  density  gradient  and  the  assumption  of  a  linear  gas  velocity  distribution.  It  can  be  shown 
that  a  zero  gas  density  gradient  implies  a  linear  gas  velocity  distribution,  but  the  reverse  is  not 
necessarily  true.  While  the  Lagrange  model  is  an  approximate  solution  to  the  Lagrange 
ballistic  problem,  it  must  be  noted  that  this  model  has  been  routinely  applied  in  the  simulation 
of  'real'  guns  in  which  the  propellant  combusts  at  a  finite  rate  during  the  ballistic  process,  in 
this  case,  it  is  assumed  that  the  unbumt  propellant  is  uniformly  distributed  throughout  the  gas, 
moves  at  the  local  gas  velocity,  and  bums  at  a  rate  dependent  on  the  space  mean  pressure. 


The  Lagrange  gradient  model  is  generally  applicable  in  cases  for  which  the  charge-to- 
mass  ratio  is  'small.'  However,  reasonable  results  are  obtained  for  C/M  in  excess  of  2.0. 
Using  the  velocity  and  pressure  distributions  of  the  Lagrange  model,  relationships  between 
breech,  space  mean,  and  base  pressures,  as  well  as  gas  and  projectile  kinetic  energies,  are 
obtained.  These  are 


where 


KEgas 


KEpRoj 


(3) 

(4) 
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is  the  ratio  of  propellant  mass  to  projectile  mass  ana  P\s  the  space  mean  pressure  of  the  gas. 
The  gradients  of  pressure  and  velocity  in  the  Lagrange  model  are  not  dependent  on  specific 
gas  properties,  other  than  C/M,  and  the  numerical  factors  which  arise  in  the  model  (1/2  and 
1/3)  are  constants  independent  of  system  parameters. 

4.  PIDDUCK-KENT  GRADIENT  MODEL 

The  Pidduck-Kent  model  is  a  special  solution  of  the  continuity  and  momentum  equations 
for  the  gaseous  propelling  charge  in  the  Lagrange  gun.  This  special  solution  was  first 
published  by  Love  and  Pidduck  (1S22).  A  more  detailed  presentation  was  latter  published  by 
Kent  (1 936)  and,  subsequently,  an  extended  treatment  was  published  by  Vinti  and  Kravitz 
(1949).  The  Pidduck-Kent  model  addresses  the  nonuniformity  of  gas  density  which  is 
neglected  in  the  Lagrange  model  and  which  becomes  increasingly  important  with  increasing 
C/M.  This  is  accomplished  by  assuming  that  each  element  of  the  gas  evolves  adiabaticaliy 
from  its  initial  state  with  all  gas  elements  following  the  same  adiabatic.  In  the  resulting 
solution,  it  is  found  that,  in  its  initial  state,  the  gas  density,  pressure,  and  temperature  are 
decreasing  functions  of  distance  along  the  tube  from  the  breech  to  the  projectile  base.  A 
detailed  derivation  of  the  Pidduck-Kent  solution  is  presented  in  Appendix  B.  Pressure  and 
energy  relationships  analogous  to  Equations  1-3  are  obtained  in  the  Pidduck-Kent  gradient 
model.  These  are 

P BREECH  *  PBASe{  1  "  aO  "  *  1  *  (5) 

p  *  paAsef 1  +  4 

l  5 

Kf:a*fKEpR(Ji)  (7) 

where  n  is  the  polytropic  index, 

1 

(7-1)’ 
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and  a 0  and  5  are  parameters  arising  in  the  Pidduck-Kent  nlution  which  are  related  by 


1_ 

5 


1 

2/7  +  3 


2(n  +  1) 
£ 


(8) 


In  the  limit  as  e  (i.e.,  C/M)  approaches  zero,  8  approaches  3,  and  a0  approaches 

£ 

T2uT7TjJ 


(see  Appendix  B)  such  that  Equations  5-7  reduce  to  the  Lagrange  model,  Equations  1-3. 

In  contrast  with  the  Lagrange  model,  the  parameters  of  the  Pidduck-Kent  solution,  a^  and 
5,  are  functions  of  C/M  and  the  ratio  of  specific  heats.  Similarly,  the  distributions  of 
thermodynamic  quantities  are  dependent  on  both  C/M  and  the  ratio  of  specific  heats.  The 
calculation  of  ag  and  5  for  a  specific  ballistic  system  has  been  simplified  by  the  development 
of  a  least  squares  fit  of  computed  values  of  5  using  a  polynomial  function  in  e  and  y  (Grollman 
and  Baer  1970).  (The  least  squares  fit  of  5  was  completed  for  1.2  £  y  5  1.3,  however,  values 
of  8  for  y  =  1.4  are  within  1-2%  of  values  obtained  using  the  tables  of  Vinti  and  Kravitz  [1949], 
see  Appendix  D).  Equation  8  is  then  used  to  calculate 

It  has  been  generally  supposed,  though  never  proven,  that  the  Pidduck-Kent  solution  of 
the  Lagrange  problem  approaches  the  wave  dynamic  solution  for  sufficiently  large  projectile 
travel.  Additional  discussion  of  this  issue  can  be  found  in  Appendices  B  and  E. 

5.  PROBLEM  DESCRIPTION 

The  Lagrange  ballistic  problem  was  chosen  as  the  framework  for  evaluation  of  the  range 
of  applicability  of  the  Lagrange  and  Pidduck-Kent  gradient  models.  This  physically  simple 
configuration  was  chosen  to  facilitate  direct  comparison  of  lumped  parameter  solutions  based 
on  these  gradient  models  with  one-dimensional  wave  dynamic  solutions.  Not  only  is  the 
Lagrange  problem  the  framework  for  the  development  of  these  gradient  models,  but,  by 
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assuming  a  prebumed  propellant,  differences  in  lumped  parameter  and  wave  dynamic 
solutions  arising  from  physical  models  of  component  processes  (i.e.,  ignition,  combustion, 
drag,  etc.)  are  avoided. 

In  order  to  investigate  the  effects  of  gas  molecular  weight,  three  propelling  gas 
compositions  were  selected.  These  are  hydrogen  gas,  which  has  been  heated  by  an  external 
source;  the  combustion  products  of  a  mixture  of  100%  hydrogen  peroxide  and  octane  with  an 
oxidizer-to-fue!  ratio  of  approximately  10  to  1;  and  the  combustion  products  of  JA2  solid 
propellant.  The  thermochemical  properties  of  these  gases,  calculated  using  the  Blake  code 
(Freedman  1982)  are  given  in  Table  1.  The  initial  gas  pressure  was  chosen  to  be  500  MPa  in 
all  cases.  The  initial  temperature  of  the  hydrogen  gas  was  chosen  to  be  the  flame 
temperature  of  JA2,  while  the  flame  temperature  of  the  octane-hydrogen  peroxide  mixture  is 
somewhat  lower. 

The  constant  bore  diameter  of  the  Lagrange  gun  was  chosen  to  be  120  mm.  The 
projectile  mass  is  9.0  kg.  Calculations  were  performed  for  the  C/M  of  0.5,  1.0,  3.0,  and  8.0 
with  the  gas  properties  presented  in  Table  1.  This  leads  to  different  initial  chamber  volumes 
in  each  case  as  indicated  in  Table  2.  Shot  start  pressure,  bore  resistance,  and  gas  pressure 
ahead  of  the  projectile  were  assumed  to  be  zero  and  heat  transfer  was  neglected  in  all 
calculations. 

The  XNOVAKTC  code  (Gough  1990)  was  used  for  one-dimensional  simulations  of  the 
Lagrange  gun.  while  the  IBHVG2  code  (Anderson  and  Fickie  1987)  was  used  in 
thermodynamic  simulations  with  both  the  Lagrange  and  Pidduck-Kent  (Grollman  and  Baer 
1970)  gradient  mode!  options  of  this  latter  code.  It  is  noted  that  the  analytic  solution  to  the 
Lagrange  problem,  as  presented  in  Appendices  A  (Lagrange  Gradient)  and  B  (Pidduck-Kent 
Gradient),  could  have  been  used  in  place  of  IBHVG2.  However,  this  latter  computer  model 
has  been  extensively  tested  and  is  readily  accessible. 
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Table  1.  Propelling  Gas  Properties 


Hydrogen 

Octane/ 

Peroxide 

JA2 

Initial  Pressure  (MPa) 

500 

500 

500 

Initial  Temperature  (K) 

3,410 

3,017 

3,410 

Initial  Density  (g/cm3) 

3.1204  x  10'2 

2.4892  x  10’1 

3.0465  x  10'1 

Ratio  of  Specific  Heats 

1.4 

1.2348 

1.225 

Covolume  (cm3/g) 

3.695 

1.044 

0.996 

Molecular  Weight  (g/mole) 

2.0 

16.873 

24.8 

Impetus  (J/g) 

14,175.4 

1,486.7 

1,140 

Specific  Energy  (J/g) 

35,438.5 

6,331.8 

5,066.7 

Gas  Sound  Speed  (m/s) 

5,036 

1,831 

1,699 

Ideal  Gas  Sound  Speed  (m/s) 

4,454 

1,355 

1,183 

Table  2.  Chamber  Volumes 


C/M  Mass  of  Gas 

(><g) 

Chamber  Volume  (liters) 

Octane/ 

Hydrogen  Peroxide  JA2 

0.5  4.5 

1.0  9.0 

3.0  27.0 

8.0  72.0 

144.2  18.08  14.74 

288.4  36.16  29.48 

865.2  108.48  88.44 

2,307.2  289.28  235.84 
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6.  RESULTS  AND  DISCUSSION 


The  results  of  the  interior  ballistic  simulations  are  presented  as  graphs  of  dimensionless 
projectile  velocity  (1 TTp )  vs.  dimensionless  projection  travel  (Tip,  see  Equations  E-21  and  E-22) 
for  C/M  of  0.5,  1.0,  3.0,  and  8.O.*  The  trajectories  for  the  hydrogen  propelling  gas  are 
presented  in  Figures  3  and  4.  At  a  C/M  of  0.5,  the  Lagrange  and  Pidduck-Kent  trajectories 
coincide  with  the  trajectory  obtained  from  the  one-dimensional  simulation.  At  C/M  of  1 .0,  the 
one-dimensional  and  Pidduck-Kent  trajectories  are  nearly  indistinguishable,  however,  the 
Lagrange  trajectory  has  begun  to  diverge.  At  C/M  of  3.0,  the  Lagrange  trajectory  has  clearly 
diverged  from  the  one-dimensional  trajectory  by  about  5%  for  dimensionless  projectile  travels 
greater  than  about  15,  which  corresponds  to  an  expansion  ratio, 


of  7.19  (see  Table  3).  From  Equation  3,  we  see  that,  for  C/M  =  3.0,  the  Lagrange  model 
predicts  that  the  gas  and  projectile  kinetic  energies  are  equal.  In  contrast,  the  Pidduck-Kent 
model  (Equation  7)  predicts  that  the  gas  kinetic  energy  is  about  85%  of  the  projectile  kinetic 
energy  (see  Table  4).  This  means  that  the  Lagrange  model  over-predicts  the  gas  kinetic 
energy  for  large  C/M,  resulting  in  an  under-prediction  of  projectile  velocity. 

The  Pidduck-Kent  trajectory  initially  diverges  from  the  one-dimensional  trajectory  for 
dimensionless  projectile  travels  between  1  and  10  (C/M  =  3.0);  however,  the  two  trajectories 
then  converge  for  xp>20.  This  behavior  is  related  to  the  wave  dynamics  in  the  gas  during  the 
early  portion  of  the  ballistic  process  (i.e.,  low  expansion  ratios).  The  energy  partition  for  the 
one-dimensional  solution,  for  C/M  =  3.0,  is  presented  In  Figure  5  as  a  function  of 
dimensionless  projectile  travel.  During  the  early  travel,  tha  wave  dynamics  of  the  system  are 
reflected  in  the  gas  kinetic  energy.  However,  the  key  feature  of  the  energy  partition  is  the 
relationship  between  gas  and  projectile  kinetic  energies.  Up  to  a  dimensionless  travel  of 

'The  introduction  of  dimensionless  variables  facilitates  the  comparison  of  results  (or  propelling  gases  o( 
substantially  different  sound  speeds  and  molecular  weights.  The  specific  choice  of  dimensionless  variables  is 
suggested  by  Siegel  (1979)  based  on  a  characteristic  solution  of  the  preburnt  propellant  ideal  gas  gun. 
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Hydrogen 


Figure  3.  Projectile  Trajectories  for  Hydrogen  Gas;  C/M  =  0.5,  1.0,  and  3.0. 


Hyd/og*n 


Figure  4.  Proiertiie  Trajectories  for  Hydrogen  Gas;  C/M  »  0.5,  1.0,  and  8.0. 
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Hydrogen 

cm=3 


Figure  5.  Energy  Partition  for  Hydrogen  Gas;  C/M  =  3.0. 


Table  3.  Relationship  Between  Expansion  Ratio  and  Dimensionless  Projectile  Travel 
for  Hydrogen 


_ 

C/M  =  0.5 

Expansion  Ratio 

C/M  =  1.0  C/M  =  3.0 

C/M  =  8.0 

5 

13.39 

7.19 

3.06 

1.77 

10 

25.77 

13.39 

5.13 

2.55 

15 

38.16 

19.58 

7.19 

3.32 

20 

50.54 

25.77 

9.26 

4.10 

25 

62.93 

31.96 

11.32 

4.87 

30 

75.32 

38.16 

13.39 

5.64 

40 

100.09 

50.54 

17.51 

7.19 

50 

124.86 

62.93 

21.64 

8.74 

100 

248.72 

124.86 

42.49 

16.48 
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Table  4.  Pidduck-Kent  and  Lagrar^e  Parameters;  y  ■-  1 .4 


C/M 

5 

■ 

kegas/kepfoj 

C/5M  C/3M  1  -  D 

PBREECk/PBASE 
(1  -  ao)-(n+1)  1  +  C/2M 

1  -  D 

0.5 

3.1272 

0.0604 

0.1599 

0.1667 

0.1591 

1 .2436 

1.25 

1 .2543 

1.0 

3.2315 

0.1055 

0.3095 

0.3333 

0.3137 

1.4775 

1.50 

1.4848 

3.0 

3.5288 

0.2174 

0.8502 

1.0000 

0.8281 

2.3581 

2.50 

2.3412 

8.0 

3.9482 

0.3447 

2.0262 

2.b667 

2.0705 

4.3895 

5.00 

4.3156 

Note:  Pidduck-Kent  parameters,  S  and  a^.  are  obtained  from  IBt  1VG2  calculations  based  on  Grollman  and  Baer 
(1970)  least-squares  fit. 

about  10,  the  gas  and  projectile  kinetic  energies  are  approximately  equal.  Above  10,  they 
diverge  and  approach  a  nearly  constant  ratio  of  gas  to  projectile  kinetic  energy  of 
approximately  0.83,  which  is  reflected  in  the  projectile  trajectory  (see  Figure  3).  Similar  results 
are  observed  for  the  C/M  =  8.0  case  in  Figures  4  and  6.  However,  the  gas  kinetic  energy  is 
greater  than  the  projectile  kinetic  energy  in  this  case  and,  due  to  the  substantially  increased 
initial  chamber  length,  a  dimensionless  travel  of  30  is  not  quite  in  the  asymptotic  (isentropic) 
region  of  the  projectile  trajectory.  However,  from  Figure  6,  it  can  be  seen  that  the  ratio  of  gas 
to  projectile  k»oetic  energy  approaches  a  constant  value,  for  large  dimensionless  travel,  which 
is  close  to  the  value  obtained  with  the  Pidduck-Kent  model  (see  Table  4).  We  also  note  that, 
despite  the  early  behavior  of  the  Lagrange  trajectory  (i.e.,  the  Lagrange  and  one-dimensional 
trajectories  are  quite  close  for  dimensionless  travels  up  to  30),  the  Lagrange  trajectory 
diverges  from  the  one-dimensional  solution  for  larger  expansion  ratios  as  the  Pidduck-Kent 
and  one-dimensional  solutions  converge. 

The  results  for  JA2  and  the  octane-hydrogen  peroxide  mixture  are  presented  in 
Figures  7-10.  In  general,  these  results  are  very  similar  to  those  already  discussed  for  the 
hydrogen  case  when  presented  in  terms  of  dimensionless  projectile  travel  and  velocity.  All 
trajectories  are  nearly  identical  for  C/M  of  0.5  and  1.0.  For  C/M  =  3.0,  the  Pidduck-Kent 
trajectories  diverge  above  the  one-dimensional  trajectories  during  early  travel  and  converge 
for  xp>25,  while  for  C/M  =  8.0  and  x^>30  ,  the  Pidduck-Kent  solutions  have  not  reached  the 
asymptotic  region  and  lie  above  the  one-dimensional  trajectories.  The  Lagrange  trajectories 
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JA2 


Figure  8.  Projectile  Trajectories  for  JA2  Combustion  Products;  C/M  =  0.5,  1. 


Octane/Hydrogen  Perioxide 


Figure  9.  Projectile  Trajectories  for  Octane/Hydroqen  Peroxide  Combustion 
C/M  =  0.5.  1.0,  and  3.0. 


0,  and  8.0. 


Products; 
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Octan* Hydrogen  ParoxWa 


Figure  10.  Projectile  Trajectories  for  Qctane/Hydroqen  Peroxide  Combustion  Products; 

C/M  =  0,5.  1.0,  and  8.0. 

for  C'M  of  3.0  and  8.0  lie  below  the  one-dimensional  trajectories  and,  in  the  C/M  =  8.0  case, 
the  Lagrange  solutions  are  expected  to  diverge  with  increasing  travel. 

7.  CONCLUSIONS 

The  applicability  of  the  Lagrange  and  Pidduck-Kent  gradient  models  in  the  simulation  of 
the  interior  ballistic  process  for  gases  of  various  molecular  weights  has  been  explored  using 
the  Lagrange  ballistic  problem  as  the  basic  framework  for  the  investigation.  The  results 
indicate  that  gas  molecular  weight  and  sound  speed  do  not  impact  the  validity  of  these 
gradient  models,  as  was  anticipated  from  consideration  of  the  theoretical  development  of 
these  models.  All  computed  trajectories  for  C/M  of  0.5  and  1 .0  are  nearly  identical.  The 
Lagrange  model  appears  to  provide  reasonably  accurate  results  up  to  C/M  of  about  3.0  {e.g., 
velocities  within  about  5%  of  the  one-dimensional  solution  at  C/M  =  3.0  with  increasing 
accuracy  for  lower  C/M).  Using  the  Pidduck-Kent  model,  the  calculated  projectile  velocity 
approaches  the  one-dimensional  solution  with  increasing  projectile  travel  for  C/M  of  3.0  and 
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1 


8.0.  The  results  also  indicate  that  the  Pidduck-Kent  and  one-dimensional  solutions  converge 
for  expansion  ratios  characteristic  of  real  guns. 
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APPENDIX  A: 

THE  LAGRANGE  GRADIENT  MODEL 
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Intentionally  left  blank. 
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The  Lagrange  gradient  mode!  was  developed  as  a  solution  to  the  Lagrange  ballistic 
problem,  which  is  described  in  the  text  and  in  Appendix  8.  This  gradient  model  is  widely  used 
in  interior  ballistic  simulations  of  real  guns,  for  which  the  original  assumption  of  a  prebumed 
propellant  gas  must  be  relaxed.  While  the  assumptions  used  in  the  derivation  of  the  Lagrange 
gradient  model  are  valid  only  for  low  charge-to-mass  ratio  (C/M),  it  has  been  found  in  practical 
applications  that  the  model  provides  reasonable  accuracy  for  C/M  as  high  as  3.0  (e.g., 
velocities  within  about  5%  of  a  one-dimensional  solution  at  C/M  =  3.0  with  increasing  accuracy 
for  lower  C/M). 

A  derivation  of  the  Lagrange  gradient  model  is  presented  here  primarily  for  completeness. 
However,  concepts  used  in  the  derivation  of  the  Lagrange  gradient  are  also  applied  in  the 
development  of  the  Pidduck-Kent  solution  in  Appendix  B.  Therefore,  inclusion  of  the 
derivation  provides  an  easy  reference  for  the  reader.  This  particular  derivation  of  the 
Lagrange  gradient  model  was  previously  published  by  Morrison  and  Coffee  (1990). 


The  equations  of  motion  governing  the  motion  of  the  gas  in  a  constant  area  gun  in  region 
from  breech  to  projectile  base  are 


3p 

9f 


=  0  , 


(A-1) 


3v  +  v  3v  dP 

dt  dx  dx 


(A-2) 


with  the  boundary  conditions 


x BREECH ~  0  > 
XBASE  ~  y 
v  BREECH  =  0  - 


(A-3) 

(A-4) 

(A-5) 
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(A-6) 


'BASE 


=  ^  =  un 


dt 


where  up  is  the  velocity  of  the  projectile. 

In  the  development  of  the  Lagrange  pressure  gradient  "it  is  assumed  ail  the  propellant 
charge  (C)  is  in  gaseous  form  at  the  time  considered."  However,  "the  theory  applies  without 
alteration  if  it  is  assumed  that,"  prior  to  consumption  of  aJI  the  propellant  charge,  "the  unburnt 
charge  moves  with  the  gas,  the  distribution  of  the  solid  along  the  bore  being  the  same  as  the 
distribution  of  gas  (Comer  1950).*  The  gas  is  assumed  to  be  inviscid  and  heat  loss  to  the 
walls  is  neglected  such  that  the  fir  w  is  isentropic. 

We  now  assume  that  the  density  of  the  gas  (or  gas  plus  unbumt  charge)  is  uniform  over 
the  region  behind  the  projectile,  i.e., 


dp 

dx 


a  0  . 


(A-7) 


We  then  obtain  from  Eq.  A-1, 


or 


0 


3v  1  9p 

Hx  pH  ’ 


(A-8) 


Assuming  a  constant  bore  area,  Ag,  (i.e.,  no  chambrage)  and  noting  that  p  =  C/Agy, 
Eq.  A-8  becomes 


dv  1  dy  up 

__  n  _ _  .  m 

dx  y  dt  y 


(A-9) 
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Integrating  Eq.  A-9  over  the  region  [0,*],  we  obtain 


v(x) 


u„ 


(A- 10) 


Corner  notes  that  the  term  “Lagrange  approximation"  is  applied  to  Eq.  A-7  or  A-10  and  the 
Eq.  A-7  leads  to  Equation  A-10,  but  it  is  not  true  that  Eq.  A-10  necessarily  implies  Eq.  A-7 
(Comer  1950). 

Substituting  Eq.  A-10  in  Eq.  A-2,  we  have 


where  up  a  du^/dt,  or 


v.'  ) 


K In 


i  a  p 

pJx 


_ .  (  \ 
dP  m  _  Cup  x 

dX  Ag 


(A-11) 


Integrating  on  [x,y],  and  noting  that 


(pbase  -  pres)  A  a 


P( 0)  =  P BREECH  ’ 


and 


p(y)  -  pbase  • 

where  PRE5  is  the  bore  resistance  pressure  and  M  is  the  projectile  mass,  we  obtain 
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p(x)  51  PBASE  +  2jj(PBASE  ~  Pf>Es) 


(A- 12) 


For  x  *  0. 


P BREECH  '  PS/SE  +  ( PBA5E  ~  PRES ) 


(A- 13) 


The  space  mean  pressure  is  defined  by 


P  •  1  f* P{x)dx  , 


(A- 14) 


and,  upon  substituting  Eq.  A-12,  we  obtain 


p  ”  pa4SE  +  ym^Pbase  ~  PR£s) 


(A-15) 


The  kinetic  energy  associated  with  the  motion  of  the  gas  is 


kegas  -  J^ylAspv2dx 


(A-16) 


Using  Eq.  A- 10,  we  obtain 


kegas  *  \  Asp 

2  r 


lx3 

3 


or 


K^gas  ■  —  Cu 


& 


(A- 17) 


26 


APPENDIX  3: 


THE  PIDDUCK-KENT  SPECIAL  SOLUTION  FOR  THE  MOTION 
OF  THE  POWDER  GAS  IN  A  GUN 


Intentionally  left  blank. 
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The  Pidduck-Kent  gradient  mode!,  like  the  Lagrange  model,  is  an  approximate  solution  of 
the  continuity  and  momentum  equations  for  the  gaseous  propelling  charge  in  the  Lagrange 
ballistic  problem.  This  solution  was  first  suggested  by  Love  and  Pidduck  (1922).  Kent  (1936) 
noted  an  error  in  Pidduck’s  solution  and  presented  solutions  for  both  ideai  and  nonideal 
gases.  Vinti  and  Kravitz  (1949)  have  published  a  detailed  treatment  of  the  Pidduck-Kent 
gradient  model,  along  with  tables  for  use  in  computing  the  parameters  which  arise  in  the 
model.  Corner  (1950)  also  discusses  the  development  of  the  Pidduck-Kent  gradient  model 

The  Lagrange  ballistic  problem  is  an  idealization  of  the  interior  ballistic  process  in  a  gun. 

It  is  assumed  that  the  gun  is  a  right  circular  cylinder  (i.e.,  the  chamber  and  bore  diameter  are 
equal)  closed  at  one  end  by  the  breech.  It  is  also  assumed  that  prior  to  the  start  of  projectile 
motion,  the  propellant  is  instantaneously  burned,  such  that  the  chamber  is  initially  filled  with  a 
quiescent  gas  of  uniform  pressure,  density,  and  temperature.’  The  gas  is  assumed  to  be 
inviscid  and  heat  loss  to  the  walls  is  neglected,  such  that  the  flow  is  isentropic. 

At  time  f  =  0,  the  initial  pressure  and  temperature  distributions  (as  in  the  case  of  the 
Lagrange  gradient  model),  as  well  as  the  initial  density  distribution  of  the  Pidduck-Kent 
solution  are  nonuniform,  decreasing  monotonicaliy  from  the  breech  to  the  base  of  the 
projectile.  Fixed  gradients  (i.e.,  distributions)  of  system  properties  throughout  the  interior 
ballistic  process  are  an  inherent  characteristic  of  gradient  models.  These  nonequilibrium  initial 
distributions  are  physically  unrealistic,  but  the  spatial  averages  of  these  distributions  must  be 
consistent  with  the  initial  equilibrium  conditions  of  the  Lagrange  problem.  However,  as 
demonstrated  computationally  in  this  work,  for  finite  charge-to-mass  ratios,  the  wave  dynamic 
solution  appears  to  approach  the  special  solution  for  large  projectile  travel,  as  has  been 
traditionally  assumed. 

In  this  Appendix,  a  detailed  development  of  the  Pidduck-Kent  gradient  model  is  presented. 
This  derivation  follows  closely  that  of  Vinti  and  Kravitz  (1949).  However,  the  notation  is 
modified  slightly  for  clarity. 


Although  the  Lagrange  and  Pidduck-Kent  gradient  models  arise  in  the  solution  of  an  idealized  problem,  they  are 
routinely  applied  in  the  simulation  of  real  gun  systems.  In  such  cases,  the  gas  is  replaced  by  a  “fluid"  consisting 
of  a  mixture  of  the  combustion  gas  and  burning  solid  propellant  grains.  The  density  of  this  fluid  is  based  on  the 
combined  masses  of  the  gas  and  solid  phases. 
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A  Lagrangian  formulation  is  used  in  the  development  of  the  Pidduck-Kent  gradient  model. 
The  gas  column  between  the  breech  and  the  projectile  base  is  initially  divided  along  the  axial 
direction  into  a  very  large  number  of  thin,  disk-shaped  elements.  A  disk  may  expand  (or 
contract)  axially,  but  the  mass  of  gas  in  each  disk  is  constant  (i.e.,  there  is  no  flow  of  gas 
between  adjacent  disks).  (Eventually,  the  number  of  disks  will  be  permitted  to  become  infinite, 
[i.e.,  the  thickness  of  the  disks  will  be  allowed  to  approach  zero]  such  that  the  solution  is 
continuous.)  The  thermodynamic  quantities  for  the  gas  (i.e.,  pressure,  density,  etc.)  are 
assumed  to  be  uniform  within  each  element,  but  are  assumed  to  vary  from  element  to  element 
in  the  axial  direction.  Let  x  be  defined  as  the  distance  of  the  midpoint  of  a  given  gas  element 
from  the  breech  and  let  y  be  the  position  of  the  projectile  base  along  the  barrel  with  respect  to 
the  breech  at  some  arbitrary  time,  t,  as  shown  in  Figure  B-1.  We  define 

x(t  =  0)*xQ  (B-1) 

where  the  subscript  "0“  denotes  the  value  of  a  given  quantity  at  t  =  0.  The  gas  is  assumed  to 
obey  the  Nobel-Able  equation  of  state,  such  that  for  a  gas  element  initially  located  at  x0 


Po(* o.O) 


1 

Po(*o*0) 


RT0{x0, 0)  , 


(B-2) 


where  P0(x0,0)  is  the  pressure  of  the  gas  in  the  element  located  at  x0  at  t  =  0,  p0(x,0)  is  the 
gas  density,  and  r0(x0,0)  is  the  gas  temperature.  As  is  implied  by  Eq.  B-2,  the  gas  covolume, 
specific  heats,  impetus,  etc.  are  assumed  to  be  constants  for  all  elements  from  the  breech  to 
the  projectile  base.  As  in  the  development  of  the  Lagrange  gradient  model,  space  mean 
thermodynamic  quantities  are  defined  as  spatial  averages  on  [0,y]  and,  for  consistency,  the 
initial  values  of  these  space  mean  quantities  must  correspond  to  the  uniform  initial  properties 
of  the  gas  in  the  Lagrange  problem.  This  will  be  addressed  in  more  detail  later. 

As  the  projectile  accelerates  down  the  barrel,  each  gas  element  is  assumed  to  expand 
isentropically, 


30 


Time  t  =  0.0 


where  the  adiabatic  constant  K(xQ)  is  a  function  of  the  initial  position  of  the  given  gas  element. 
In  order  to  make  the  problem  more  tractable,  it  is  assumed  that  all  gas  elements  follow  the 
same  adiabatic,  K,  such  that 


Consider  now  an  element  of  gas  initially  located  at  xQ  with  thickness  Ax0  such  that  the 
mass  of  the  gas  in  the  element  is  p(x0,0)  As  Ax0,  where  AB  is  the  cross-sectional  area  of  the 
bore.  At  some  later  time,  f,  the  gas  element  will  have  moved  to  a  position  x(x0,  /)  and  will 
have  expanded  such  that  its  thickness  is  A x[xQ,t}  as  shown  in  Figure  B-2.  Taking  the  limit  as 
A x{xQ,t)  approaches  zero,  the  continuity  equation  for  this  gas  element  is  then 

po(xo,0  )dx0  =  p  (xQ,t)dx(x0,t)  (B-5) 

and  the  equation  of  motion  is 
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Time  t 


Figure  B-2.  Evolution  of  a  Gas  Element  From  Time  f  =  0  to  Time  t. 


P  (x0,t) 


d2x(x0,t)  BP(x0,t ) 


dX 


(B-6) 


which  are  the  continuity  and  momentum  equations  in  a  lagrangian  frame  of  reference. 


Following  Vinti  and  Kravitz  (1949),  we  introduce  a  new  variable 


r(x0.f)  ■  x{x0,t)  -  t]  j*°p0(x£  ,0)dx£  , 


(B-7) 


such  that 


dz(xQlt) 


dx(xQ,t) 

8*0 


-  riPo(xO'°)  . 


(B-8) 


and 


z(x0,t)  _  dx(x0lt) 

dt  =  Tt 


(B-9) 
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Rewriting  the  continuity  equation  (Eq.  B-5) 


dx(x0,t)  _  po(*o-°) 
3*o  p(*o-0 


(B-10) 


and  substituting  in  Eq.  B-8,  we  have 


Bz(x0,t) 

3*o 


Po(^o-O) 


1 

P(*  o.O 


(B-1 1) 


We  now  define 


W(x0it) 


1 

p(*o.O 


(B-1 2) 


such  that  Eq.  B-4  becomes 


P(x0,t)  =  K[W(x0,t)r 


(B-1 3) 


and  Eq.  B-1 1  becomes 


dz{x0,t) 

3*o 


Po(^.0)[lV(*o.O]'1  • 


(B-1 4) 


Using  Eq.  B-9,  B-10,  and  B-13,  the  equation  of  motion  can  be  rewritten  as 

B2z(x0lt)  K 

9 12  Po(xo>0)  3x0 


(B-1 5) 
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where  we  have  made  use  of 


BP{x0,t)  _  dP{x0,t)  dx0  _  p(x0,f)  dP(x0,t) 
3xq  Jx  p0(Xq,Q)  Jx~q 


(B-16) 


We  now  attempt  a  solution  of  Eq.  B-15  by  separation  of  variables  assuming  a  solution  of 
the  form 


z{x0,t)  =  f(x0)<S>(t)  . 


(B-1 7) 


At  t  =  0, 


2(x0.0)  =  ^{^0)^(0)  = 


xQ-n^Po(xi,o)dxi 


0(0) 


(B-1 8) 


such  that 


f(*o)  =  *o  -  ti  I/  Po(xo  •  0)  dx o 


(B-1 9) 


and 


0(0) =  1  . 


(B-20) 


Substituting  in  Eq.  B-1 4,  we  obtain 


dz{xQ,t) 

d*o 


df{x  0) 
dx0 


®{t)  *  Po(x0,0)[W(xQ,t)) 


-1 


and  using  Eq,  B-1 9  to  evaluate 


df{x0)  . 

- we  have 

dx  o 
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where 

a^.o,  _  .  iWj,,  ,  ,Pq(j,0i0)1  ,  (B.25) 

ix0  d/(x0)  dx0  a/(r,)  l  01  0  ’> 

such  that 

-yKlW0(x0,0)V-2  dWo(xo,0)_a 
f(*  o)  a'<*0) 
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fb a  fiYo)  -  Yo-nC/AB 


(B-31) 


where 


S(n,a0)  s  JJ  (i  -  aQiiz(x0))nd\i(x0). 


(B-40) 


It  is  convenient  to  rewrite  Eq.  B-40  at  this  point.  Let 

q(p)  =  30  \i2  (Xq) 


such  that 

^(xo)= 

2/a0q(n) 


Substituting  in  Eq.  B-40,  we  obtain 


S(n,a0)  =  foa°q(M)’2(l  -qM)ndq(n)  (B-41) 


which  gives  us 

S(n,a o)  =  —1^  B  J'n  ♦  1 ) - 2Fi f-«4 4 :ao]  (B-42) 

2yc?o  \  ^  ^  J 

1  13 

where  Bao(_,  n  +  1 )  is  the  incomplete  beta  function  and  2^i  (-n, _. ;ao )  is  Gauss’ 

hypergeometric  function.  Eq.  B-42  will  prove  useful  later  in  our  development  of  the 
Pidduck-Kent  gradient  model. 

Eq.  B-40  could  be  utilized  to  calculate  the  parameter  a^  given  a  detailed  description  of  the 
initial  state  of  the  system.  However,  Vinti  and  Kravitz  (1949)  have  developed  a  more  general 
equation  from  which  a q  can  be  determined.  Note  that  at  the  base  of  the  projectile,  the 
equation  of  motion  becomes 
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which,  with  the  use  of  Eq.  B-27  and  B-33,  becomes 

y 

=  Agfb  lVo(0,0)[1  -  a0)~  .  (B-46) 

Multiplying  both  sides  of  Eq.  B-46  by  S{n,a 0)  and  using  Eq.  B-39,  we  obtain 

.  i 

_il^a0(1  -  a0)  I*  S(n,a 0)  =  (B-47)  - 

which  can  be  used  to  calculate  the  parameter  a0  given  only  the  ratio  of  specific  heats  of  the 
gas  and  the  propelling  charge  to  projectile  mass  ratio,  C/M. 

Ratio  of  Breech  to  Base  Pressure 

The  determination  of  the  ratio  o*  breech  pressure  to  projectile  bass  pressure  is  now 
straightforward.  Using  Eq.  B-27  and  B-33,  we  have 

i 

Wo(yC.0)  -  w0(0,0)[1  -  ao]'-1  (B-48) 

or,  dividing  both  sides  by  <&(t)  and  using  Eq.  B-21, 

i 

W(y0,t)  »  W(0,t)  [1  -  a0]~  .  (B-49) 

We  now  make  use  of  Eq.  B-13  to  obtain 

4r~h’  (1  -a0  (B-50) 
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which  demonstrates  that  at  all  times  during  gas  expansion  the  ratio  of  breech  to  projectile 
base  pressure  is  a  constant  dependent  on  7  and  a§. 


Space  Mean  Pressure,  Density,  and  Temperature 

As  noted  earlier,  space  mean  thermodynamic  quantities  are  defined  as  spatial  averages 
on  [0 ,yj  and,  for  consistency,  the  initial  values  of  these  quantities  must  correspond  to  the 
uniform  initial  state  of  the  gas  in  the  Lagrange  problem.  As  a  first  step,  we  write  the  equations 
for  the  spatial  distributions  of  pressure,  temperature,  and  density  in  the  gas.  From  Eq.  B-12, 
B-13,  B-21,  B-27,  B-33,  and  B-34,  we  have 


P(*o.O  -  Po(0.0)[l  -  a042(x0)]',  +  1[<X>(f)]-7  (B-52) 

and, 

T(*0,f)  -  T0 (0,0)  [  1  -  a0n2(x0}]  [O(0]1_Y  •  (B-53) 

Consider  first  the  gas  density.  We  rewrite  Eq.  B-51  in  the  form 

W{xQ,t)  -  Wo(0,0)  [l  -  aoU2^)]"  O(f)'1  .  (B-54) 

Performing  a  spatial  average  over  0  s  p,  s  1,  we  obtain 

W)  -  ^(O.OHd^f)]-1  [1  .  a0n2(*0)]V(*o) 

or,  using  Eq.  B-40, 

W)  -  W0{0.0)(  O(f)]-1  2FJ-n';3;a0).  (B-55) 

V  2  2  J 
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This  result  can  also  be  obtained  from  Eq.  B-7  and  B-17, 


’(yo*0  -  =  y{t)  -  h/>j*0',o)d*'  =  y(t)  -  Ct\ia8  (b-56) 


where  we  have  noted  that  x(y0,  0  s  MO  is  the  position  of  the  projectile  base  as  a  function  of 
time.  Since  the  average  gas  density,  p(f),  is  the  gas  mass  divided  by  the  gas  volume, 


4,0(0  =  -n  . 

Mp(0  J 


(B-57) 


where  y(t)  is  the  position  of  the  projectile  base  at  some  time,  f,  after  the  start  of  projectile 
motion.  Rearranging  and  using  Eq.  B-39  and  B-42,  we  have 


=  -n  =^o(0 
P(0  J  c 


(B-58) 


W(t)  -  WoiO.OJjOC/)]  -1  2  F i  -n .1  ;l,a0 


(B-59) 


If  we  approach  this  problem  in  a  slightly  different  manner  using  Eq.  B-17  and  B-31,  we  obtain 


zLo .0  =  4,0(0  =  (y0  -  riCMs)0{f) 


40(f)  =  .£  _L_  -  n  o(o 
Mp(0)  J 


42' 


Again  rearranging  and  using  Eq.  B-39  and  B-42,  we  have 


P(0)  )  C 


or,  using  Eq.  B-12 


mo  )mt)}-'  =  IV0(C,0)[<I><  f)] -1 2  :  )  * 


which,  with  Eq.  B-59,  gives  us 


[  w<0]-1  •[w(0)]-'*[t) 


(B-60) 


(B-61) 


Thus,  the  spatial  average  of  W(t)  which  is  a  function  of  the  gas  density,  has  the  same  time 
dependence  as  W(x0,  f)  for  each  individual  gas  element,  which  is  given  by  Eq.  B-21.  This 
result  might  well  have  been  anticipated  as  a  resuit  of  the  assumption  that  all  gas  elements 
follow  the  same  adiabatic. 

We  now  consider  the  spatial  average  of  the  gas  temperature,  which  will  be  defined  in 
terms  of  the  internal  energy  of  the  gas,  i.e., 


£/(*)■  C 


(B-62) 


where  £r,(f)  is  the  internal  energy  of  the  gas  at  any  time,  t,  and  c^is  the  specific  heat  at 
constant  volume.  Since  energy,  like  mass,  is  an  extensive  quantity,  we  have,  using  Eq.  B-29, 


Ej(t)  =  AgCy  T(x0,f)  p(x0,f)  dx{x0,t) 


(B-63) 


where  the  integration  is  now  over  the  gas  column  at  some  arbitrary  time,  t.  We  first  use 
Eq.  B-7  to  change  variables 


-(x0,f)  =  x(x0,0  -  q  Joxo  p (xg  ,0 ) dXQ  , 
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where,  from  Eq.  B-56, 


*b  s  2(y0,t)  =  y{t)  -  nC/Ab  . 


Using  Eq.  B-10  in  B-11,  we  have 


dx(x0,t)  =  [1  -  tip^o-O]"1  dz(x0,t) 


such  that  Eq.  B-63  becomes 


f/(0  -AbC*  j*b  T (x0,t) 


1 


P(*o •t) 


-  T1 


-1 


dz(x0,t) 


We  now  define  a  new  variable 


\i{XQ)  s 


z[x0,t)  _  z(x0,t)  _  f(x0)<t>(t)  ~  f(x 0) 
z{y0J)  zT~  TyJWo  fb 


which  is  identical  to  Eq.  B-34,  our  original  definition  of  We  now  note  that 


dz(x0lt)  =  zbdv.(x0)  =  44>(f)«W*0) 


where  Eq.  B-63  and  Eq.  B-65  are  taken  at  a  fixed  time,  f,  and  thus 


and 


H(0)  =  0, 

n(y0)  = 1 


as  before.  Using  Eq.  B-51,  B-53,  and  B-67  in  Eq.  B-65,  we  obtain 


(B-64) 


(B-65) 


(B-66) 


(B-67) 
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Ei(t)=ABc,To(0,0) 


1 


Po(°-0) 


-  *1 


|-i 

[i  -  a0ii2(x0)}n*' dvixo) 


or,  using  Eq.  B-17,  B-31,  and  B-40, 


Ei(t)  =  0,70(0,0) 


1 


Po  (0»0) 


-  n 


-i 


ABfb  2F,[-n  -  1.i;l;a0)[®(0]Hr 


Using  Eq.  B-39  and  B-42,  we  note  that 


ABfb 


1 


Po(O.O) 


-  T1 


1-1 


2Fi  (~n,i.;£;a0) 


such  that  Eq.  B-68  becomes 


W)  -  C^TofO.OH^f)]1^ 


2Fi(-n  -  l,l;|;a0) 
-  - 

2  ( -/i.  — ;  — ;  a0 ) 


and,  using  Eq.  B-62,  we  obtain  the  space  mean  temperature 


T(0  =  T0(0,0)  [cD(f))1^ 


2 Fj ( -n  -  1,l;|;a0) 
2  Fi  ( -n,  1 ;  £ ;  a0 ) 


or,  using  Eq.  B-53,  and  noting  that  p(0)  =  0, 


(B-68) 


(B-69) 


(B-70) 


(B-71) 
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T{/)  =  T(O.f) 


(B-72) 


zF,(-n-ui;l;ao) 
2^1  (~n'-^>^>ao) 


In  order  to  determine  the  space  mean  pressure,  we  note  that 


P(t) 


Pit) 


~  T| 


=  R  T(0 


(B-73) 


or,  using  Eq.  S-12, 


P(t)  =  R  T (t)  W(t) 


(B-74) 


Combining  Eq,  B-55  and  B-72  and  using  Eq.  B-21,  we  obtain 


1  .  3 


P(t)  =  m(0,t)  W(O.f)  2F,(-n  -  1,A;^;ao) 


or 


W)  =  P(0,t)  2Fi(-n-h±;l;a0)  . 


(B-75) 


Now  we  divide  Eq.  B-75  by  P(y0,f)  and  use  Eq.  B-50  to  obtain 


(B-76) 


From  Eq.  B-47,  we  have 

(1  -  ao)-<n+1)  = 


2 [n  +  1 )  a0  2  F)  ( -n,  ~ ;  ;  aQ ) 


(B-77) 
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which  in  Eq.  B-76  gives  us 


P(t) 
PiYo.i ) 


e 

2(n  +  1 ) 


2  p\  ( ~n  - 1  ;  -g- ;  ao ) 

2Fi  (~n'4:'§  :a°^ 


(B-78) 


It  can  be  shown  that  (see  Appendix  C) 


Fi  (  -n  -  1  J.;ao)  =  (1  -  3$) 


n*  1  .  2  3o  (n  +  1 ) 


2Fi(-n. 


3 . 5 
2  ’  2 


;  ) 


(B-79) 


or,  using  Eq.  B-77 


Pin  i  1  .  3  .a  >  2(n  + 

2  Fi  ( -n  -  1 .  — ,  — ,  So )  = - 

2d.  e 


1 )  a<3 


2  cto  (r>  +  1 ) 

3 


2  F1  ( -n,  1 ;  2 ;  a0 ) 


(B-80) 


which,  in  Eq.  B-73,  gives  us 


-fllL  =  i  +  £ 

F(y0,f)  3 


2^1  (_n>-2  :4:ao) 


(B-81) 


or 


-  i  ♦£ 

F(y0,O  5 


(B-82) 


where 


2 

8 


1  2  Fi  ( _n  ’  § ;  § ;  ^ ) 
2fi  (-n,l;£;a0) 


(B-83) 
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Noting  that 


is  real  and  finite  on  0  5  <,  1  for  n  >  1  of  practical  interest,  we  see  from  Eq.  B-47  that  =  0 

corresponds  to  e  =  0  and  that  Sq  =  1  corresponds  to  e  -»  ».  Using  the  definition  of  the 
Gauss'  hypergeometric  function  (see  Appendix  C)  we  see  from  Eq.  B-83  that  8  =  3  for 
e  =  Sq  =  0  (as  expected)  and  5  =  2n  +  3  for  -  1.  Thus,  in  the  limit  of  zero  C/M  the 
Pidduck-Kent  solution  reduces  to  that  of  Lagrange. 

As  is  also  demonstrated  in  Appendix  C,  a  relationship  can  be  derived  among  the 
parameters  a 0,  e,  n,  and  8, 

1  „  1 
8  =  2  n  +  3 

which  can  be  used  to  determine  8  for  a  given  e  and  n  once  a 0  has  been  calculated.  For  small 
e,  it  is  easily  shown  from  Eq.  B-84  that  ag  approaches  e/2(n  +  1). 

Projectile  Velocity 

The  velocity  of  an  element  of  gas  is  given  by  Eq.  B-9 

dx(x0,t)  _  dz(x0,t) 

Tt  =  Ti 


_i_ 

30 


2 (n  +  1) 


(B-84) 


and,  using  Eq.  B-17,  we  obtain 

u(x0,t) 


dx(x0,t) 

87 


(B-85) 


from  Eq.  B-23,  we  have 

•  (8-86) 
dt 2 
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We  now  define 

¥(0  ■ 

df 


such  that  Eq.  B-86  can  be  rewritten  in  the  form 

d2^{t)  _  <mt)  d$>(t)  _  lt)  <W(t)  _  flr<D(nr 
cf<D(/)  dt  d<L>(t) 


which*  upon  integrating  with  respect  to  <5(f)  and  noting  that  0(0)  =  1,  gives  us 


and  from  Eq.  B-85 


[^(f)]2  =  f^!lT  =  _iL  {i  -t‘Kor(T‘1)} 
l  dt  ) 


u2(x0,t)  «  C/(X0)]2  JLL{l  -  [0(f)]-^-i)} 


From  Eq.  B-28  and  B-33,  we  have 

b-  ■  2yKa°  4  nvoio.oir* 

(T-1  )<£ 


where,  from  Eq.  B-13, 

/<S  P0(0.0)[W0(0,0)H 

such  that 

B  =  P0(0,0)[W0(0,0)]-1 

(Y  -  1 ) ^ 


or 


2730 

S- - PT0(0,0) 

(Y*  1)4 


(8-87) 


(B-88) 


(B-89) 


(B-90) 


49 


which  gives  us 


50 


in  the  form  presented  by  Siegel  (1979).  We  now  use  Eq.  B-72  to  evaluate  T(0,0)  in  terms  of 
the  initial  energy  of  the  gas.  For  t  =  0,  we  have 


T(0)  =  T0(0,0)  - -  (B-96) 

2Fi('n’2^’>ao) 


which,  using  Eq.  B-78  and  8-80,  gives  us 


T(0)  =  To{0,0) 


2(rt  +  1}3o 

( 

1  +  - 

e 

6 

V  ) 

or 


T0(0,0)  =  T(0) 


1 

2  (n  +  1 )  a0 


(B-97) 


Substituting  Eq.  B-97  in  Eq.  B-95  and  noting  that 


RT(  0) 
7-1 


*  E,( 0)  a  Eq  , 


(B-98) 


the  initial  specific  internal  energy  of  the  gas,  and 


27 
7  -  1 


=  2  (n  +  1)  , 


we  obtain 


u2(y0,t)  - 


2£0e 

1 


1  - 

1 

(B-99) 
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which  is  solution  to  the  preburned  propellant,  nonideal  gas  gun  assuming  the  Pidduck-Kent 
gradient  model,  which  is  discussed  in  Appendix  E.  In  the  limit  of  very  large  projectile  travel, 
Eq.  B-95  becomes 


(|2  4YRTo(0,0)ao 
UP  (7-  I)2 


(B-100) 


which,  using  Eq.  B-97  and  B-98,  gives  us 


Ay  RT(0) 
(Y-1)2 


(  e  }  5 

e  +  6  ^  2  (n  +  1 ) 


(B-101) 


It  has  been  generally  supposed  that  the  Pidduck-Kent  solution  of  the  Lagrange  problem 
approaches  the  wave  dynamic  solution  at  sufficiently  large  projectile  travel.  Corner  (1950)  has 
developed  an  expression  of  the  form  of  Eq.  B-100  for  the  maximum  possible  velocity  for  the 
Pidduck-Kent  gradient  model  for  the  case  of  infinite  C/M.  Recalling  that  =  1  for  e  °°,  we 
have 

i 

umax  =  -JL.  (yRT0(0,0))2  .  (B-102) 


which  is  identical  in  form  to  the  escape  velocity  obtained  in  the  characteristic  solution  of  the 
Lagrange  problem.  However,  in  the  Lagrange  problem  the  initial  gas  pressure,  temperature 
and  density  are  uniform  throughout  the  chamber.  Thus,  in  the  initial,  physical  state  of  the  gas, 
the  temperature  at  the  breech  and  the  space  mean  temperature  are  equivalent.  In  the 
Pidduck-Kent  solution,  the  initial  pressure,  density,  and  temperature  distributions  decrease 
monotonically  from  the  breech  to  the  base  of  the  projectile,  a  physically  unrealistic  initial  state 
for  the  system.  Thus,  the  space  mean  temperature,  from  which  the  specific  internal  energy  of 
the  gas  is  determined,  is  lower  than  the  initial  breech  temperature.  Therefore,  for  equivalent 
jas  internal  specific  energies,  the  Pidduck-Kent  solution  results  in  a  maximum  (escape) 
velocity  greater  than  that  obtained  from  the  characteristic  solution  of  the  Lagrange  problem  for 
infinite  C/M  and  infinite  travel.  From  Eq.  B-101  we  obtain 


52 


Comer  (1950)  has  mistakenly  equated  the  breech  temperature  in  Eq.  8-102  with  the  initial 
temperature  of  the  Lagrange  problem  and  then  concluded  that,  at  least  for  infinite  C/M,  the 
Pidduck-Kent  and  wave  dynamic  solutions  yield  the  same  maximum  velocity.  As 
demonstrated  here,  the  Pidduck-Kent  maximum  velocity  is  actually  somewhat  higher  than  the 
escape  velocity  for  y  >  1 .  This  difference  is  about  8%  for  y  =  1 2,  14%  for  y  =  1 .4,  and  20% 
for  y  =  1 .667.  It  would  appear  that,  at  least  for  the  infinite  C/M  case,  ihe  Pidduck-Kent 
solution  does  not  approach  the  wave  dynamic  solution  at  large  projectile  travel,  as  has  been 
traditionally  assumed.  However,  as  shown  in  the  comparisons  of  calculated  trajectories  for 
finite  C/M,  the  Pidduck-Kent  solution  does  appear  to  approach  the  wave  dynamic  solution  at 
sufficiently  large  projectile  travel. 


Gas  Kinetic  Enerc 


The  kinetic  energy  of  the  propelling  gas  at  an  arbitrary  time,  t,  is  defined  by 


2 

dx(xQ,t)  . 


(B-105) 


Using  Eq.  B-64,  B-66,  B-67,  and  B-85  we  obtain 


kegas  -  4  ABfb*{t)  (I  u2(*o.t)  --■■  ■  -  it 
2  J0  p(*o.O 


d\i(x0)  . 


(B-106) 
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Evaluation  of  the  Function  f(xc): 


The  equations  developed  thus  far  are  adequate  for  solution  for  the  interior  ballistic  problem 
to  determine  projectile  velocity  as  a  function  of  travel.  In  addition,  the  breech,  projectile  base, 
and  space  mean  values  of  the  thermodynamic  functions  (i.e.,  pressure,  temperature,  etc.) 
can  be  calculated  as  functions  of  projectile  travel.  Given  the  values  of  these  functions  at 
the  breech  prior  to  the  start  of  projectile  motion,  Eq.  B-51-B-53  could  then  be  used  to 
evaluate  the  thermodynamic  functions  at  any  arbitrary  position  along  the  barrel  and,  thus,  the 
spatial  distribution  of  the  functions,  for  any  given  projectile  travel.  However,  Eq. 

8-51  -B-53  oxpress  the  thermodynamic  functions  in  terms  of  a  dimensionless  variable  ^(Xg) 
and  it  is,  therefore,  necessary  to  determine  the  relationship  between  x0  and  p(x0)  if  we  wish  to 
investigate  the  behavior  of  spatial  variation  of  these  functions  in  the  Pidduck-Kent  solution. 

The  dimensionless  function  ^(Xq)  is  defined  by  Eq.  B-34, 


where 


and 


such  that 


H(*o)  3 


f(*  o) 


f{x o)  -  *0  -  n  J0X°Po(^o/-0)c/xo 


M*o)  = 


f{yo)  -  /o  -  nC/AB 
*°  “  ‘n  /0X° p  o  ( • 0  > 


(B-34) 


(B-19) 


(8-31) 


(B-112) 


We  must  now  evaluate  the  integral,  C{Xq), 


C(x 0)  =  Ag  |0x0p(xo/,0 )dx0' 


(B-113) 
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such  that 


56 


and 


<?(  0)  =0. 


Substituting,  W8  obtain 

' 

C(xQ)  =  ABfbW0(0,0)  1%  2  J^^qbi)  2(1  -  g(n))',d<7ai)J  (B-116 


or,  using  Eq.  C-17,  we  have 

' 

-1  -1 

C(x0)  -  ABfbWQ(0,0)  \i(xQ)  l(a0n2(x0>)  2  f^*0' q(n)  2(1  -  q(\i))n  dqfc) 
or 

C(x0)  =  Ae^Wo(0,0)  fi(x0)  (B-117) 


Using  Eq.  B-117  in  Eq.  B-114  and  rearranging  terms,  we  obtain 

*o  -  nOfoJkj1  +  n2frJ-n,l;j;30n2(x0)  J  (B-ii8) 


which  is  the  desired  functional  relationship  between  x0  and  p.(xQ).  The  numerical  evaluation  of 
this  function  is  discussed  in  Appendix  D. 


Summar 


The  fundamental  assumption  of  the  Pidduck-Kent  gradient  model  is  that  all  elements  of 
gas  in  a  Lagrange  gun  expand  along  the  same  adiabatic  ( K )  such  that 


P o(*o-0) 


Po(*o«0) 


(B-2) 
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which  is  rewritten  by  defining  W0(x0,0)  such  that 


P0(xQlO)  =  KlW(x0,0)F  .  (8-13 


A  new  variable  is  introduced 


*(*0.0  =  *(*0.0  ■  ^ $0° Po (x0  *°) dxo  -  f(Xo)*(t)  (B-7),  (B-17 


and  a  solution  of  the  momentum  equation  is  attempted  by  separation  of  variables.  An 
expression  for  the  total  mass  of  propelling  gas  in  the  system  is  obtained  by  integrating  the  gas 
density  distribution  over  the  length  of  the  chamber, 


C  =  ABfbW0( 0,0)  2Fi(-n.^;*;a0) 


(B-39),  (B-42 


where  fb  =  f{y0),  n  is  the  polytropic  index, 


1 

Y-  1  * 


and  is  a  parameter  which  must  be  determined  in  order  to  satisfy  this  relation.  After  some 
manipulation,  a  more  basic  relation  is  found. 


2nao(1  -  a0)-<n+1>  2F,  (-n,-l;|;a0)  =  e  (B-47 

where  e  is  the  gas  propelling  charge  to  projectile  mass  ratio. 

It  is  shown  than  each  element  of  gas  evolves  in  time  according  to  the  relation 

[ IVfxo.f)]'1  =  [ Wq(Xo.O)]'1  3>{f)  . 


(B-21 
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The  gas  kinetic  energy  is  found  to  be 


KE, 


GAS 


fO 

UJ 


KE 


PROJ 


(B-110) 


such  that 


KE- 


TOTAL  = 


KE \ 


PROJ  ■ 


(B-111) 


which  are  easily  seen  to  reduce  to  the  Lagrange  result  in  the  limit  of  zero  C/M. 


The  projectile  velocity  as  a  function  of  expansion  ratio  is  determined  to  be 


u  (y0’t) * 


2  E0e 

1  +£ 


1  - 


( 1  V-i 

£f 

V.  '  / 


(B-99) 


which  is  the  solution  of  the  preburned  propellant,  nonideal  gas  gun  problem  as  can  be 
obtained  directly  from  energy  considerations.  Finally,  it  is  shown  that  the  maximum  projectile 
velocity  obtained  from  the  Pidduck-Kent  solution  of  the  infinite  travel,  infinite  C/M  Lagrange 
problem  is  given  by 


P-K 

Viax 


"  3y  -  1  n 

1 

1 

L  2y  J 

L  7  -  1 


(y  Fmo)) 


3y  -  l" 

L  2y  J 

J escape 


(B-103,  B- 1 04) 


where  uescapo  is  the  maximum  velocity  obtained  in  the  exact  solution  of  this  problem. 
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APPENDIX  C: 

DISCUSSION  OF  INTEGRALS  ARISING  IN  PIDDUCK-KENT 
SPECIAL  SOLUTION  OF  THE  LAGRANGE  PROBLEM 
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Intentionally  left  blank. 
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Several  integrals  arise  in  the  development  of  the  Pidduck-Kent  special  solution  of  the 
Lagrange  ballistic  problem  which  can  be  related  to  the  Gauss’  hypergeometric  functions  and 
incomplete  "beta"  functions  (Arfken  1970;  Abramovitz  and  Stegun).  These  relationships  prove 
useful  in  the  derivation  of  the  Pidduck-Kent  solution  and  in  determining  its  limiting  behavior. 
The  integrals  of  interest  are  presented  in  Appendix  B  in  Eq.  B-38,  B-68,  and  B-107, 


h  s  /J  (l  -  a0p2)/’cfp  ,*  (C-1) 

h>  =  JJ  (l  -  a0p2),?  +  1  d\i  ,  (C-2) 

and 

Is  m  (l  -  aoP2)np2cfp  .  (C-3) 

In  Eq.  C-1,  we  define  a  new  variable  q  such  that 

q  s  aoP2  (C-4) 

1 

H  -  ( q/ao )*  (C'5) 


d\i  -  _J__g  2  dq 
2yao 


(C-6) 


q(o)  =  0  (C-7) 

g(1)  =  So  (C-8) 


*  The  integral  in  Eq.  8-112  defining  C(x^  is  closely  related  to  /?  and  will  not  be  considered  separately. 
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giving  us 


30 


/*[i  -  q}nq 


(C-9) 


Similarly,  Eq.  C-2  and  C-3  become 


-q\Mq~?dq 


and 


k 


3 


1 

- qYq^dq . 


(C-10) 


(C-11) 


Gauss’  hypergeometric  function  is  defined  by  (Arfken  1970;  Abramovitz  and  Stegun) 


2  (a,b\c\z) 


r(c)  f i 
r(b)r(c  -  b)  Jo 


(C-12) 


wherb  T(c)  is  the  gamma  function.  For  c  -  b  *  1, 


and, 


such  that 


r(c  -  b)  =  r(i)  -  i 


r(0  .  r(b  +  i)  _  brjb)  _  b 
'r(b)  r(b)  r(b) 


(013) 


(014) 


2F,  (a,b;b+1  ;z)  =  b  t^(  1  -  tzya  dt . 


(C-15) 


We  now  define  a  new  variable,  q  =  tz,  such  that 

2^(a,b;b  +  1 ; z)  =  J.  -  0‘a^  (C-16) 


64 


1  -  2  M 

4  =  So  Ba0-5--/,+2  * 


(C-23) 


1  '■5  7  3  A 

B*  !,#,  +  1  • 


(C-24) 


We  now  integrate  Eq.  C-2  by  parts  to  obtain 


huVL[  1  "  So^2]^1  I  +  2(n  +  1 ) ao J^1  [ 1  -  ao p. 2 ] n 2 c/p.  ,  (C-25) 


h.  m  (  1  ~  ao) n+1  +  2 ( n  +  1 ) 3q/3  . 


(C-26) 


Alternately,  we  could  proceed  by  rewriting  the  integrand  in  Eq.  C-2  to  obtain, 


>2  "  fo  I1  "  ao^2rdu  "  ^  £<1  -  a0p.2)i7p2cTp.  ,  (C-27) 


/2  =  /i  -  a0/3  . 


(C-28) 


Equating  Eq.  C-26  and  C-23.  we  have 


h  -ao*3  ={1  _ao)n+1  +  2  ( n  +  1 )  a0/3  , 


which  can  be  rewritten  as 


1  ^  (i-aor1 

2n  +  3  Sq  a0 


(C-29) 


(C-30) 


From  Eq.  B-47,  we  have 


0  -  So)' 


rt+1  _  2(n  +  1)30/! 


(C-31) 
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and  substituting  in  Eq.  C-30,  we  obtain 


or 


„  1  1  _  2(/J  +  1) 

/}  2n  +  3  2q  e 


(C-32) 


1  2fi 

3  2^l(^4;|;ao) 


1 

2n  +  3 


2(/>  +  1) 
e 


(C-33) 


which  is  a  relationship  between  two  contiguous  hypergeometric  functions  (Arfken  1970; 
Abramovitz  and  Stegun).  Using  Eq.  B-81  to  identify  the  left-hand  side  of  Eq.  C-33  as  1,  we 
have 


1  =  1  1  _  2(n  +  1) 

8  2n  +  3  3q  e 


(C-34) 


Eq.  C-28.  also  an  equation  for  contiguous  functions,  can  be  rewritten  as 


2  F\  ( ~n  - 1  •  ;  2o ) 

2  ( -n,  A  ;  £ ;  a0 } 


a0  2  (  n  •  i  a0 ) 

3  aF1{-/».^;|;^)) 


fo 

T 


(C-35) 


Following  the  arguments  of  Vinti  and  Kravitz  (1949),  we  first  note  that,  from  Eq.  B-50, 


P{y0’t) 

P(o.t) 


[1  -  aop*1) 


(C-36) 
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it  is  dear  that  0<  1.  Since  n  s - and  y  >  1,  then  n  is  a  real,  positive  number. 

y  -  1 

Thus,  if  the  pressure  ratio  given  in  Eq.  B-50  is  to  be  a  real  number  and  P[o, f)  £  P(yQ,f),  it  must 
follow  that  0  £  Sq  5  1.  We  also  note  that  the  integrals  of  interest,  Eq.s  C-1,  C-2,  and  C-3  are 
real  and  finite  for  n  >  o  and  0<  aQ<  1.  Then,  from  Eq.  B-47 


e  -  2  (n*  1 )  a0  ( 1  -  2F,{-n,^-a0) 


(C-37) 


we  see  that  »  0  corresponds  to  the  limit  of  zero  C/M  while  3q  =  1  corresponds  to  the  limit 
of  infinite  C/M.  Using  Eq.  C-1 5  to  evaluate  the  hypergeometric  function  for  Sq  =  0,  we  have 


2 F]  ( a,b\b  *  1  ;0)  =  bfj  t ^  dt=  tb\  - 1 


(C-38) 


For  3g  =  1,  the  integrals  of  interest  reduce  to  regular  (complete)  ’beta"  functions,  giving  us 


(C-39) 


(C-40) 
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Recalling  that 


1  (3  ^  ,rir(n*1’ 

4  /  ■  5\ 

^  >  r  n  + 1 
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rUp) 


(  g  V  i  .  , 

2  n  +  _  r,  +  _  r  /7+_ 

l  2  A  2  ]  t  2 


1  2Fi(-/7,|;|;^)  /3 

3  'i 


(C-41 ) 


(C-42) 


we  obtain  for  Sq  =  0 


2  _  2 
5  =  3 


(C-43) 


and  for  Sq  =  1 , 


B  _ ,  n  +  1 

1  _  l2 

1  "71 - 

Bf^l 


2  r?  +  3 


(C-44) 


This  latter  relation  can  be  obtained  directly  from  Eq.  C-34.  However,  the  zero  C/M  limit  is 
shown  to  be  a  direct  result  of  the  solution,  in  contrast  to  the  argument  used  by  Vinti  and 
Kravitz  (1949)  involving  expansion  of  e  in  Eq.  C-31in  terms  of  ^  and  substitution  of  the  result 
into  Eq.  C-34. 

The  ratio  Iji ,  arises  in  Eq.  B-72  relating  the  space  mean  gas  temperature  to  the 
temperature  at  the  breech.  For  =  1,  we  obtain 


;2(ao  a  T) 

/i(ao  *  1) 


(n  +  l)  _  2{n  +  l)  ^  2y 

'  3l  2  n  +  3  7  -  1 

n  +  _ 

2 


(C  ') 


which  is  the  factor  arising  in  Eq.  B-101  and  B-102  relating  the  maximum  velocity  from  the 
Pidduck-Kent  solution  for  the  infinite  travel,  infinite  C/M  Lagrange  gun  problem  and  the  escape 
velocity  obtained  from  the  exact  solution  of  this  problem. 
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APPENDIX  D: 


NUMERICAL  EVALUATION  OF  THE  PIDDUCK-KENT  PARAMETERS 


Intentionally  left  blank. 
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The  evaluation  of  the  Pidduck-Kent  parameters,  %  and  8,  defined  by  Eq.  B-47  (or  C-37) 
and  Eq.  B-82, 

e  -  2(n  +  1)ao(1  -  ao)-(n+1)  (D-l) 


and 


where 


1  1  1  2(n  +  1) 

S’  2n  +  3  ao  e 


(D-2) 


2Fi(-/7,^:|:a°)  =  JJfl  -  ao ji2)n dp. 


(D-3) 


has  been  considered  previously  by  Vinti  and  Kravitz  (1949)  and  by  Grollman  and  Baer  (1970). 
Vinti  and  Kravitz  developed  an  algebraic  expression  for  8  as  a  function  of  n, 


1  _  1 
S’  2/7  +  3 


1  +  a(e) n 


'l  +  C1(e,/7)p(e)n> 

1  +  C-i  (e.n)  n 
\  / 


(D-4) 


where  the  coefficients  a(e),  P(eJ,  and  C^e.n)  are  themselves  functions  of  e  and  n.  These 
coefficients  are  tabulated  (Vinti  and  Kravitz  1949)  for  selected  values  of  e  and  n  such  that  the 
value  of  the  coefficients,  and  thus,  8  can  be  determined  for  arbitrary  e  (0  £  e  ^  10)  and 
n  ( 0.5  <  n<  5  [i.e.,  1.2<,  y<,  3\)  by  interpolation.  While  the  computed  values  of  8  obtained 
using  the  tables  of  Vinti  and  Kravitz  are  quite  accurate,  this  technique  is  cumbersome  and  not 
well-suited  for  application  in  computer  simulations.  In  contrast,  Grollman  and  Baer  developed 
extensive  tables  for  5  as  a  function  of  e  and  n,  for  0  <  e  £  20,  and  1.2  £  y  S  1.3,  by  numerical 
integration  to  obtain  e  as  a  function  of  for  given  n  (see  Eq.  D-1  and  D-3)  and  use  of 
Eq.  D-2  to  obtain  8.  A  least  squares  fit  of  the  tabulated  values  of  8  was  then  generated  using 
a  fifth  order  polynomial  function  of  e  and  y, 


5 


3.0 


5 

+  z 
/»1 


«/■ + 


(D-5) 
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This  expression  is  easily  applied  in  interior  ballistic  simulations,  and  has  served  as  the  primary 
means  for  implementation  of  the  Pidduck-Kent  gradient  model.  However,  in  using  Eq.  D-5  to 
compute  8,  errors  on  the  order  of  a  few  percent  are  introduced.  Additionally,  the  range  of 
applicability  of  this  equation  is  limited,  particularly  in  view  of  the  properties  of  some  novel 
propellants  being  considered  in  conjunction  with  the  electrothermal-chemical  gun. 

Our  objective,  then,  is  to  develop  a  general  technique  for  calculation  of  the  Pidduck-Kent 
parameters,  a 0  and  5,  which  is  valid  for  arbitrary  e  and  5,  and  which  ii?  readily  applicable  in 
computer  simulations.  As  a  by-product,  we  also  develop  techniques  to  evaluate  the  other 
functions  encountered  in  the  Pidduck-Kent  solution  of  the  Lagrange  problem. 

The  approach  selected  to  develop  a  generally  applicable  method  for  determination  of  the 
Pidduck-Kent  parameters  is  direct  solution  of  Eq.  D-1  for  a§.  The  hypergeometric  function 
may  be  expressed  as  an  infinite  series  of  the  form,  i 


9  2Fi 


T(k  +  <$>)T(k  •+  y) 


r(<t>)r(¥)  o 


where  the  series  converges  for  0  -  y  -  4>  >  0  and  |x|  <  1.  For  9  =  y  +  1,  the  Gauss 
hypergeometric  function  reduces  to  the  incomplete  "beta"  function  for  0  <  x  <,  1  and  the 
regular  (complete)  "beta"  function  for  x  =  1.  For  the  functions  of  interest  in  the  Pidduck-Kent 
solution,  0  =  y  +  1  and  <j>  =  -n  or  -  (n  +  1)  where  9,  y  and  m  are  positive  numbers  such  that 
0  —  y  —  (fi  >  0.  We  also  note  that  x  =  a0  with  0  <  <  1,  such  that  the  series  in  Eq.  D-6 

converges  for  the  values  of  parameters  considered  here. 

We  are  generally  interested  in  Gauss  hypergeometric  functions  of  the  form 


2^(1  -  b,a\a  +  1  ;x) 


T(a  +  1)  ^  r(k  ♦  1  -  b)f(k  +  a)  xk 

r(1  -  b)  F(a)  o  r(k  +  1  +  a)  k\ 


(D-7) 


We  will  now  simplify  this  expression  using  the  properties  of  the  "gamma"  function.  First,  we 
note  that 


r(a  +  1)  ar(a)  0 
r(a)  "  r(a) 


and 


T(k  +  a)  _  r(/c  +  a)  1 

T(k  +  1  +  a)  (k  +  a)T(A  +  a)  k  +  a 


We  now  rewrite  Eq.  D-7,  to  obtain 


2^(1  -  b,a;a  ♦  1;x)  «  1  +  a  X  _ _  (D-10) 

k. i  T(1  -  b)  (*  +  a)/t!  v  ' 


The  ratio  of  "gamma"  functions  can  be  expanded  to  give  us 
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such  that 
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(D-12) 


Eq.  D-1  can  now  be  solved  iteratively  using  the  Newton-Raphson  technique.  We  define 


Giao.me)  =  2(n  +  1)a0(1  -  ao)-(n+1)  2^i  -  e  -  0  (D-13) 
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such  that 


where 


and 


dG{aQ,n,£) 
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An  initial  estimate  for  ag  is  obtained  using  the  Groliman  and  Baer  (1970)  formula 

(Eq.  D-5)  to  calculate  8(1)  and  Eq.  D-2  to  calculate  a^1*.  Successive  values  of  e$  are  then 

calculated  iteratively  using  the  Newton-Raphson  equation 
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until  the  desired  accuracy  is  achieved. 

In  order  to  improve  convergence  of  the  infinite  series  in  Eq.  D-7  for  0.65  <xs  1.0,  we 
consider  the  integral  representation  of  the  hypergeometric  function 

2F,{1  -  b,a\a  +  1;x)  =  (1  -  0b_1  dt .  (D-18) 
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Changing  variables,  we  define 
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which  (Abramovitz  and  Stegun)  we  might  have  used  as  the  starting  point  to  obtain  the  desired 
relation,  Eq.  D-19.  From  the  symmetry  of  the  'beta"  function,  i.e., 


B(a,b)  -  *  B (b,a)  , 

r(a  +  o) 


we  see  that  Eq.  D-21  can  be  rewritten  as 


B x(a,b)  =  B[a,b)  -  B Ux(b,a) 


and  that  B (b,a)  could  be  similarly  replaced  in  Eq.  D-19. 


Returning  now  to  Eq.  D-1S  for  the  arguments  of  interest,  we  have 


2  a0  2 


(  n  1  -3-a 
-n,  _ ,ao 


=  B 


(1 


_,n  +  1 
2 


(1  3q) 


n*  1 


n  +  1 


2  V 


1 


;  n  +  1 ;  n  +  2 ;  1 


where 


Bj-l,n  +  1 

l2 


(D-22) 


(D-23) 


-  ao  |  (D-24) 


(D-25) 


(D-26) 


78 


and 
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We  also  require  the  derivative  with  respect  to  a$. 
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Therefore,  for  3q  >  0.65,  Eq.  D-13  and  D-14  become 
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In  practice,  is  generally  less  than  0.5  for  the  values  of  e  and  y  of  interest.  However,  we  do 
utilize  this  alternate  formulation  in  comparison  of  a 0  and  8  calculated  using  the  other 
techniques  (Vinti  and  Kravitz  1949;  Grollman  and  Baer  1970). 

The  computational  technique  described  above  has  been  implemented  in  a  computer 
program,  which  is  presented  in  Appendix  F.  This  program  has  been  used  to  calculate  and 
8  for  the  values  of  e  and  n  from  the  tables  of  Vinti  and  Kravitz  (1949).  The  results  of  these 
calculations  are  presented  in  Tables  D-1-D-4,  along  with  corresponding  values  of  ag  and  8 
obtained  using  the  techniques  developed  by  Vinti  and  Kravitz  (1949)  and  Grollman  and  Baer 
(1970).  Comparing  the  values  of  2q  and  5  calculated  using  the  techniques  presented  here 
with  those  of  Vinti  and  Kravitz,  we  see  that  they  generally  agree  within  1  part  in  1 05.  We  note 
that  Vinti  and  Kravitz  (1949)  have  estimated  that  values  of  3q  and  8  calculated  from  their 
tables  have  an  accuracy  of  1  part  in  6,700  and  1  part  in  5,000,  respectively.  The 
convergence  criteria  for  the  Newton-Raphson  iteration  used  in  our  calculations  provides  an 
accuracy  of  at  least  1  part  in  105  in  the  determination  of  3q.  Therefore,  the  technique 
presented  here  is  at  least  as  accurate  as  that  of  Vinti  and  Kravitz  over  the  range  of  e  and  y 
considered.  As  noted,  this  technique  has  been  implemented  in  a  computer  program, 
therefore,  our  initial  objective  has  been  accomplished.  We  also  note  that  the  accuracy  of  the 
Grollman  and  Baer  (1970)  formula  is  somewhat  less  than  that  of  the  other  two  approaches, 
varying  fromm  about  1  part  in  300  for  charge-to-mass  ratios  around  3,  to  about  1  part  in  1000 
at  higher  and  lower  charge-to-mass  ratios. 

Grollman  and  Baer  (1970)  have  published  tables  of  the  Piddi»~:;-Kent  parameters  (6)  for 
0  S  e  £  26  in  increments  of  0.01  and  for  1.2  <  y<  1.3  in  increme.  ,iS  of  0.01.  These  results  are 
based  on  a  numerical  integration  to  obtain  e  as  a  function  of  a^  for  given  n  and  use  of 
Eq.  D-2  to  calculate  8.  The  values  of  5  at  the  tabulated  points  we.  e  obtained  by  interpolation. 
Selected  values  of  8  from  these  tables  are  presented  in  Tables  D-5-D-7,  along  with 
corresponding  values  of  8  calculated  using  the  least  squares  fit  of  Grollman  and  Baer  (1970) 
and  the  technique  developed  in  this  work  in  order  to  assess  the  accuracy  of  the  current 
computational  method  for  10  £  e  S  20.  As  can  be  seen  from  these  tables,  our  values  agree 
quite  well  with  the  tabulated  values  of  Grollman  and  Baer.  The  accuracy  of  the  Grollman  arid 
Basr  least  squares  fit  is  comparable  to  that  noted  earlier  at  lower  charge-to-mass  ratios. 

Finally,  we  note  that  the  accuracy  of  the  tabulated  values  of  5  from  Grollman  and  Baer  is 


5  parts  in  105  (based  on  the  convergence  criteria  for  their  numerical  integrator)  which  is 
consistent  with  comparisons  to  values  of  8  obtained  using  the  technique  developed  in  this 
work. 

As  an  additional  note,  we  point  out  other  applications  of  the  techniques  developed  in  this 
Appendix.  In  order  to  calculate  %  as  a  function  of  p  for  0  <  p  £  1  from  Eq.  B-1 18,  or  to 
evaluate  the  other  hypergeometric  functions  arising  in  the  Pidduck-Kent  solution  (i.e.,  Eq.  C-19 
and  C-20)  we  can  utilize  Eq.  D-12  and  D-19. 

The  computer  program  used  to  calculate  a 0  and  5  with  the  numerical  method  described 
above  has  been  incorporated  into  the  lumped  parameter,  interior  ballistic  code,  IBHVG2. 
IBHVG2  is  a  standard,  well-documented,  IB  mode!  for  solid  propellant  guns  containing  the 
Grollman  and  Baer  least  squares  fit  for  e.  The  IB  code  was  modified  to  permit  the  use  of 
either  the  least  squares  fit  or  the  numerical  technique  described  in  this  report  to  determine  ag 
and  8. 

Since  the  accuracy  of  the  least  squares  fit  decreases  with  increasing  e,  a  sample  problem 
of  interest  with  charge-to-mass  ratio  of  3.003  was  examined.  A  120-mm  gun  with  7  perf,  JA2 
propellant  was  chosen  as  a  baseline.  The  values  of  ag,  8,  maximum  breech  pressure  and 
muzzle  velocity  for  each  case  is  shown  in  Table  0-1.  Three  cases  are  presented:  (1)  least 
squares  fit  with  an  optimized  web;  (2)  numeric  method  using  the  web  determined  in  (1);  and 
(3)  numeric  method  with  an  optimized  web. 

A  comparison  of  the  three  cases  in  Table  0-8  shows  that  the  numeric  method  gives 
approximately  a  0.1%  difference  in  maximum  breech  pressure  and  a  0.03%  difference  in 
muzzle  velocity  compared  to  the  least  squares  fit  for  calculating  8  for  a  given,  optimized  web. 
However,  if  the  web  is  re-optimized,  the  differences  drop  to  0.007%  in  maximum  b:eech 
pressure  and  0.004%  in  muzzle  velocity.  Thus,  for  high  charge-to-mass  ratios  in  a  regime  of 
interest  (i.e.,  -  3  in  this  application)  differences  in  predicted  gun  performance  using  fitted  or 
numerically  calculated  values  for  8  and  ag  in  the  Pidduck-Kent  gradient  are  considered 
insignificant. 
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Table  D-1.  Comparison  of  Pidduck-Kent  Parameters  for  y=  1.2 


Vinti  and 

Gro  liman  and 

Morrison  et  al. 

Kravitz  (1949) 

Baer  (1970) 

E 

*0 

6 

% 

5 

*0 

8 

0.2 

0.01557 

3.06355 

0.01556 

3.05887 

0.01557 

3.06352 

0.4 

0.02927 

3.12164 

0.02926 

3.11479 

0.02927 

3.12159 

0.6 

0.04150 

3.17532 

0.04149 

3.16790 

0.04150 

3.17517 

0.8 

0.05255 

3.22502 

0.05252 

3.21837 

0.05255 

3.22497 

1.0 

0.06260 

3.27150 

0.06258 

3.26633 

0.06260 

3.27156 

2.0 

0.10259 

3.46897 

0.10264 

3.47325 

0.10259 

3.46887 

3.0 

0.13134 

3.62597 

0.13201 

3.63614 

0.13184 

3.62604 

4.0 

0.15481 

3.75762 

0.15502 

3.76739 

0.15481 

3.75754 

5.0 

0.17365 

3.87054 

0.17381 

3.87664 

0.17366 

3.87099 

0.18963 

3.97123 

0.18962 

3.97115 

0.18962 

3.97099 

1 

0.20343 

4.06084 

0.20328 

4.05609 

0.20342 

4.06053 

8.0 

0.21558 

4.14194 

0.21533 

4.13486 

0.21557 

4.14168 

0.22641 

4.21608 

0.22616 

4.20S39 

0.22640 

4.21593 

10.0 

0.23616 

4.28408 

0.23601 

4.28047 

0.23617 

4.28440 
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Table  D-2.  Comparison  of  Pibduck-Kent  Parameters  for  1.25 


Vinti  and 

Gro  liman  and 

Morrison  et  al. 

Kravitz  (1949) 

Baer (1970) 

£ 

*0 

6 

3d 

5 

3) 

5 

0.2 

0.01866 

3.06099 

0.01866 

3.05648 

0.01866 

3.06097 

0.4 

0.03505 

3.11672 

0.03504 

3.11009 

0.03505 

3.11666 

0.6 

0.04966 

3.16816 

0.04964 

3.16101 

0.04966 

3.16801 

0.8 

0.06281 

3.21577 

0.06279 

3.20936 

0.06281 

3.21572 

1.0 

0.07477 

3.26027 

0.07474 

3.25529 

0.07477 

3.26033 

2.0 

0.12211 

3.44893 

0.12217 

3.45312 

0.12211 

3.44884 

3.0 

0.15649 

3.59850 

0.15670 

3.60841 

0.15649 

3.59856 

4.0 

0.18335 

3.72355 

0.18360 

3.73312 

0.18334 

3.72347 

5.0 

0.20527 

3.83054 

0.20546 

3.83662 

0.20528 

3.83093 

6.0 

0.22378 

3.92579 

0.22379 

3.92591 

0.22377 

3.92553 

7.0 

0.23972 

4.01030 

0.23955 

4.00601 

0.23971 

4.01001 

C.O 

0.25370 

4.08667 

0.25342 

4.08019 

0.25369 

4.08644 

9.0 

0.26612 

4.15634 

0.26585 

4.15033 

0.26612 

4.15624 

10.0 

0.27727 

4.22015 

0.27713 
—  _ _ 

4.21718 

0.27729 

^  .22051 
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Table  D-3.  Comparison  of  Pidduck-Kent  Parameters  for  y  =  1.3333 


Vinti  and 

Grollman  and 

i  Morrison  et  al. 

Kravitz  (1949) 

Baer  (1970) 

e 

*0 

8 

^0 

5 

8 

0.2 

0.02329 

3.05716 

0.02328 

3.05288 

0.02328 

3.05713 

0.4 

0.04368 

3.10933 

0.04366 

3.10306 

0.04368 

3.10928 

0.6 

0.06179 

3.15745 

0.06176 

3.15066 

0.06179 

3.15732 

0.8 

0.07806 

3.20194 

0.07802 

3.19584 

0.07805 

3.20190 

1.0 

0.09281 

3.24349 

0.09277 

3.23872 

0.09280 

3.24354 

2.0 

0.15078 

3.41909 

0.15083 

3.42293 

0.15077 

3.41902 

3.G 

0.19244 

3.55765 

0.19267 

3.56G82 

0.19243 

3.55772 

4.0 

0.22471 

3.67305 

0.22498 

3.68175 

0.22469 

3.67295 

5.0 

0.25085 

3.77136 

0.25104 

3.77661 

0.25085 

3.77177 

6.0 

0.27279 

3.85354 

0.27275 

3.85808 

0.27276 

3.85837 

7.0 

0.29158 

3.93577 

0.29132 

3.93091 

0.29155 

3.93552 

8.0 

0.30797 

4.00526 

0.30758 

3.99821 

0.30794 

4.00510 

9.0 

0.32248 

4.06864 

0.32207 

4.06176 

0.32245 

4.06849 

10.0 

0.33545 

4.12640 

0.33518 

4.12228 

0.33544 

4.12670 
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Table  D-4.  Comparison  of  Pidduck-Kent  Parameters  for  y  =  1.5 


Vinti  and 

Kravitz  (1949) 

Groilman  and 

Baer  (1970) 

Morrison  et  al. 

e 

*0 

b 

ao 

5 

ao 

5 

0.2 

0.03097 

3.05078 

0.03096 

3.04690 

0.03097 

3.05075 

0.4 

0.05794 

3.09705 

0.05792 

3.09132 

0.05794 

3.09700 

0.6 

0.08177 

3.13964 

0.08174 

3.13341 

0.08177 

3.13954 

0.8 

0.10307 

3.17900 

0.10303 

3.17329 

0.10307 

3.17895 

1.0 

0.12230 

3.21565 

0.12225 

3.21108 

0.12230 

3.21569 

2.0 

0.19696 

3.36984 

0.19702 

3.37259 

0.19696 

3.36978 

3.0 

0.24966 

3.49055 

0.24991 

3.49745 

0.24966 

3.49060 

4.0 

0.28988 

3.59030 

0.29014 

3.59607 

0.28988 

3.59025 

5.0 

0.32208 

3.67490 

0.32217 

3.67655 

0.32209 

3.67514 

6.0 

0.34879 

3.74927 

0.34853 

3.74497 

0.34878 

3.74912 

7.0 

0.37147 

3.81492 

0.37085 

3.80567 

0.37145 

3.81466 

8.0 

0.39107 

3.87361 

0.39020 

3.86150 

0.39106 

3.87348 

9.0 

0.40829 

3.92687 

0.40732 

3.91406 

0.40828 

3.92682 

10.0 

0.42358 

3.97536 

0.42267 

3.96401 

0.42360 

3.97560 
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Table  D-5.  Comparisons  of  Pidduck-Kent  Parameter  5; 
Y  =  1 .2  and  y  -  1 .22  for  e  >  1 0. 


Y  =  12 

Y=  1.22 

Grollman  and  Baer  (1970) 

Morrison  et  al. 

Grollman  and  Baer  (1970) 

Morrison  et  a). 

E 

CALC 

FIT 

CALC 

FIT 

10 

4.2844 

4.2805 

4.28440 

4.2581 

4.2545 

4.25808 

12 

4.4073 

4.4115 

4.40727 

4.3777 

4.3824 

4.37773 

14 

4.5153 

4.5233 

4.51525 

4.4828 

4.4912 

4.48276 

16 

4.6116 

4.6128 

4.61164 

4.5764 

4.5781 

4.57642 

18 

4.69S7 

4.6836 

4.69871 

4.6610 

4.6516 

4.66095 

20 

4.7781 

4.7811 

4.77814 

4.7380 

4.7416 

4.73799 

Table  D-6.  Comparisons  of  Pidduck-Kent  Parameter  5; 
Y  =  1 .24  and  y  =  1 .26  for  e  >  1 0. 


Y  *  1.24 

Y  =  1.26 

Grollman  and  Baer  (1970) 

Morrison  et  al. 

Grollman  and  Baer  (1970) 

Morrison  et  al. 

E 

CALC 

FIT 

CALC 

FIT 

10 

4.2328 

4.2294 

4.23279 

4.2085 

4.2051 

4.20846 

12 

4.3474 

4.3541 

4.34935 

4.3221 

4.3268 

4.32208 

14 

4.4516 

4.4602 

4.45157 

4.4216 

4.4301 

4.42161 

16 

4.5426 

4.5446 

4.54263 

4.5102 

4.5121 

4.51020 

18 

4.6248 

4.6158 

4.62475 

4.5900 

4.5812 

4.59001 

20 

4.6995 

4.7035 

4.69952 

4.6626 

4.6665 

4.66264 
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Table  D-7.  Comparisons  of  Pidduck-Kent  Parameter  5; 
Y  =  1 .28  and  y  -  1 .30  for  e  >  1 0. 


Y  =  1.28 

Y  1.30  | 

Grollman  and  Baer  (1970) 

Morrison  et  al. 

Grollman  and  Baer  (1970) 

Morrison  et  al. 

e 

CALC 

FIT 

CALC 

FIT 

10 

4.1850 

4.1816 

4.18504 

4.1625 

4.1588 

4.16248 

12 

4.2958 

4.3003 

4.29584 

4.2706 

4.2746 

4.27059 

14 

4.3928 

4.4009 

4.39282 

4.3651 

4.3727 

4.36512 

16 

4.4790 

4.4606 

4.47904 

4.4491 

4.4502 

4.44909 

18 

4.5567 

4.5476  ' 

4.55666 

4.5246 

4.5151 

4.52462 

20 

4.6272 

4.6307 

4.62724 

4.5932 

4.5960 

4.59324 

Table  D-8.  Comparison  of  Influence  of  Calculating  5  and  a0  Using  Least  Squares  Fit  vs. 
Numeric  Method  on  Gun  Performance 


Calculation 

8 

20 

Max  Breech 
Pressure 
(MPa) 

Muzzle 

Velocity 

(m/s) 

Fit 

(optimized  web) 

3.62239456 

0.14475450 

690.0536 

2461.12 

Numerical  method 
(same  web) 

3.61239195 

0.14456433 

689.2609 

2460.30 

Numerical  method 
(optimized  web) 

3.61239195 

0.14456433 

690.1069 

2461.03 
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Intentionally  left  blank. 
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APPENDIX  E: 

LIMITING  VELOCITY  FOR  PREBURNED  PROPELLANT, 
NONIDEAL  GAS  (PPNIG)  GUN 


89 


The  statement  of  the  prebumed  propellant,  nonideal  gas  (PPNIG)  gun  problem  is  identical 
to  that  of  the  Lagrange  ballistic  problem.  It  is  assumed  that  the  gun  is  a  right  circular  cylinder 
(i.e.,  no  chambrage)  closed  at  one  end  by  the  breech.  (Siegel  [1979]  relaxes  this  assumption 
and  treats  preburned  propellant  gun  cases  with  and  without  chambrage.  For  our  purposes  we 
consider  only  the  latter.)  Prior  to  the  start  of  projectile  motion,  the  propellant  is  assumed  to  be 
instantaneously  burned,  such  that  the  chamber  is  initially  filled  with  a  quiescent  gas  of  uniform 
pressure,  density,  and  temperature.  The  gas  is  assumed  to  be  inviscid  and  heat  loss  to  the 
walls  is  neglected,  as  are  losses  due  to  resistive  forces  (i.e.,  the  flow  is  assumed  to  be 
isentropic).  The  solution  of  this  problem  can  be  accomplished  either  by  use  of  a  gradient 
model,  as  in  Appendices  A  and  B,  or  by  numerical  integration  of  the  equations  of  motion  for 
the  system.  If,  however,  the  objective  is  to  determine  a  limiting  velocity  for  this  problem,  it  is 
assumed  that  the  projectile  travel  and  charge-to-mass  ratio  (C/M)  are  permitted  to  approach 
infinity.  The  solution  of  the  resulting  problem  is  then  straightforward. 

Characteristic  Solution 


The  characteristic  equations  for  an  effectively  infinite  travel,  infinite  C/M  prebumed 
propellant  gun  is  presented  by  Siegel  (1979)  for  the  ideal  gas  case.  An  extension  to  the 
nonideal  (Nobel-Abel)  gas  equation  of  state  is  presented,  following  the  development  of  Siegel. 

When  the  projectile  begins  to  move,  a  rarefaction  wave  moves  back  into  the  quiescent  gas 
at  the  speed  of  sound,  c0.  If  the  chamber  is  sufficiently  long,  either  there  will  be  no  wave 
reflections  from  the  breech  or  the  reflected  wave  will  not  reach  the  projectile  base  prior  to 
muzzle  exit.  The  equation  of  state  of  the  gas  is 


(  ^ 


RTn 


the  gas  internal  specific  energy  is 


(E-1) 


(E-2) 
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and  the  gas  sound  speed  is 


c 2  _  Y  Pp  m  Y*T0 

Po(1"T1P0)  (1  -TlPo)2 


(E-3) 


We  note  that,  consistent  with  our  assumptions  regarding  the  initial  state  of  the  system,  the 
thermodynamic  properties  of  the  gas  are  uniform  throughout  the  gas  such  that  space  mean 
and  local  values  of  any  quantity  are  equal.  Since  the  gas  is  assumed  to  expand  isentropically 
as  a  result  of  projectile  motion, 


P(x.t) 


(  1 

,P  (x.t) 


=  K=P0 


r, 


The  Riemann  function,  is  defined  by 


(E-4) 


(E-5) 


where  s  denotes  constant  entropy.  Since  we  have  assumed  the  expansion  process  is 
isentropic,  this  subscript  will  be  dropped.  We  note  that  the  definition  of  the  sound  speed  has 
been  used  in  Eq.  E-5  to  relate  dP  and  dp,  i.e., 


c 


2  = 


J?)s' 


(E-6) 


Using  Eq.  E-1,  E-3,  and  E-4,  it  can  be  shown  that 


cp  =  yP0[yRTo) 


1 

1 


■\  y+1 

y 


(E-7) 
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Substituting  Eq.  E-7  in  Eq.  E-5  and  integrating,  we  obtain 


4  = 


(ynr0) 2 
yp0 


111  Jr*  i\ 

P02y  jp  [  2r  I  dP 


(E-8) 


or 


(yPTc)1 

yp0 


(p  1 

r  0 

7*1 

27 

'  2yP  " 

-s 

f  2  VP 

[P) 

Ly-1  j 

u- 1  Jcp 

(E-9) 


From  Eq.  E-3  we  note  that 


such  that 


*  c(l  -  ^P)  -  (y^t)2 


7  2  ' 


iyP t)  2  . 


(E-10) 


(E-11) 


The  identical  result  is  obtained  in  the  case  of  an  ideal  gas,  where  c2  =  tRT. 


The  characteristic  solution  for  the  system  is 


«  +  5  « 


(E-12) 


where  u  (P&  Q  =.  0.  It  is  convenient  to  evaluate  §  and  using  Eq.  E-9  which  can  be 
rewritten  in  the  form 


\  = 


'  2  1 
U  - 1 J 


1 

(Y^0}2 


(E-13) 
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such  that  Eq.  E-12  becomes 


or 


u  = 


(Y*7-0) 


f  ') 

Y-1 

1  - 

p 

27 

Pn 

- 

<  °) 

- 

(E-14) 


P 

Po 


(y  -  1  )u 

i 


2T 

T-1 


2(yRT0)2 


(E-15) 


When  the  gas  pressure  at  the  base  of  the  projectile  drops  to  2ero,  the  gas  can  do  no  more 
work  and  the  projectile  has  attained  the  maximum  possible  velocity  or  the  "escape  velocity," 

(Y RT0)*  .  (E-16) 


-“escape 


Y  "  1 


While  the  escape  velocity  does  not  represent  a  physically  attainable  velocity,  it  does  serve  as 
an  upper  velocity  limit  for  gas-driven  guns. 

The  relationship  between  projectile  velocity  and  projectile  travel  is  obtained  from  the 
equation  of  motion, 

dun 

M—P-PbAB.  (E-17) 


Evaluating  Eq.  E-15  at  the  projectile  base  and  substituting  in  Eq.  E-17,  we  have 


du„ 

M _ £  =  PnA 

dt 


o*B 


1  - 


<Y  -  1)wp 

i 

2{yRT0)1 


il 

Y-1 
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or 


Mup =  PnA 
p  dxp 


B 


(y  - 1  )>JP 

i 


il 

7-1 


2(iRT0)? 


Integrating,  we  have 


(xp  P°Ab 

Jo  M 


1 


(Y-  1  )u' 
1 

2(Y^T0)2 


27 

7-1 


1 

and  defining  a  new  variable,  z  s  up/(yR70)  2  ,  Eq  E-19  becomes 


P0A&xp 

MyRT0 


rz 

Jo 


i 

1 

1 

i - 

N 

L  l  2  J 

!  J 

7-1 


z'  dz‘ 


which  gives  us  the  desired  result 


Po^Bxp 

MyR70 


2 

Y  +  1 


2  7  +  1 


7-1  7-1 


1  - 


Y  -  1  UP 
2  1 
(Y^T0)^ 


7-1  Up 


(Y*T0) 


7+1 

7-1 


>  • 


(E-18) 


(E-19) 


(E-20) 


95 


We  now  introduce  reduced  variables  based  on  those  suggested  by  Siegel  (1979) 


— — 

Xp  =  M-(RT0 


and 

i 

u~p  =  up(yRT0)  2 

such  that  Eq.  E-20  becomes 


2 

y  +  i 


2  _  y  +  1 

y  - 1  y  - 1  L 


i  - 


i  -izlm 


-,Y  +  1 


+  1 


r  - 1 


(E-21) 


(E-22) 


(E-23) 


We  note  now  that  the  reduced  projectile  travel  can  be  rewritten  in  terms  of  the  projectile  C/M, 
e,  for  finite  chamber  length  xQ, 


P  Y(1  -  5P0) 


(E-24) 


which  makes  it  possible  to  relate  reduced  travel  to  expansion  ratio  as  defined  in 
Appendix  B 

c  _  xp  +  x0{1  '  HPo) 
f  ~><o^  -  npo) 


(E-25) 


Substituting  Eq.  E-24  in  Eq.  E-25,  we  obtain, 


ef 


(i  -  np0)*c 


(E-26) 
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We  also  note  that  the  reduced  projectile  velocity  (Eq.  E-22)  can  be  rewritten  In  terms  of  the 
initial,  specific  internal  energy  of  the  gas  ( E0 )  to  give  us 


Up  -  up  [y(y  -  1)E0]  2  . 


(E-27) 


These  equations  are  used  in  the  text  in  plots  cf  ballistic  trajectories. 


Solution  Usinq  Piciduck-Kent  Gradient  Model 


The  Pidduck-Kent  solution  to  the  Lagrange  ballistic  problem  (i.e.,  the  projectile  velocity  as 
a  function  of  gas  expansion  ratio)  is  obtained  in  Appendix  B  in  conjunction  with  the  derivation 
of  the  gradient  model.  We  show  in  this  section  that  this  solution  can  also  be  obtained  in  a 
straightforward  manner  from  the  requirement  that  energy  be  conserved.  The  initial,  total  gas 
internal  energy  is 


P  _  CRT0 

Co  —  - 


(E-28) 


ttie  initial  space  mean  pressure  is 


7T  1  J  OT 

'O  a  -  "1  R'o  *  -T - 7T- 

p o  AB*o  ~  Cti 


(E-29) 


and  the  space  mean  pressure  at  any  projectile  travel  { x ^  is 


p  ab(x0  ♦  Xp)  -  Cn 


(E-30) 


where  the  internal  energy  at  xp  is  obtained  by  noting  that,  in  the  absence  of  dissipative 
mechanisms, 

£(*p)  =  £0  -  KEPROj  -  KEgas  . 


(E-31) 
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As  demonstrated  in  Appendix  B,  for  the  Pidduck-Kent  gradient  model, 


p(^)=  1  +lr^> 


(E-32) 


KEqas  =  j'KEp 


(E-33) 


Combining  these  equations,  we  have 

^(*p)  =  ~  -p  1  +  £ 

*  V  <V 


(E-34) 


Pbitp) 


7  -  1 

'  2E0 

2 

i  +  £ 

.  s 

Ab(xo  +  •*□) 


(E-35) 


The  equation  of  motion  for  the  projectile  is  given  by  Eq.  E-17 


dun 

M — £  =  PbAB  . 
dt  b  B 


(E-17) 


Using  Eq.  C-33  in  Eq.  C-17  and  rearranging  terms,  we  have 


up^p  -  Xi_L 

P  Cfxp  2 


+  xo{  ^  “  *1  Po) 


(E-36) 


ui  du!> 

i  ♦  £ 

5 


f'p  dxp 

Jo  7  ~ - - 

XP  +  xo('  -  np0) 


(E-37) 
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from  which  we  obtain  the  desired  result 


1  +  £  L 


^  V 


(E-38) 


where  E0  is  the  initial,  specific  internal  energy  of  the  gas,  and 

c  =  xp  +  x 0 (1  -  hp0) 

£/  *0{1  -  hPo) 


(E-39) 


is  the  expansion  ratio  of  the  gas.  In  the  iimit  of  infinite  C/M  and  infinite  expansion  ratio 
(Eq.  E-38)  yields  the  maximum  projectile  velocity  for  the  Pidduck-Kent  solution 


(E-40) 


or,  recalling  that  n  =  .  J. .  , 

Y-1 


P-K 


(3y  -  1 

l  2? 


(E-41) 


Comparing  with  Eq.  E-16,  we  see  that 


P-K 


f3y  -  1 

l 


u, 


escape 


(E-42) 


as  obtained  directly  in  Appendix  B  (see  Equation  B-104)  in  the  derivation  of  the  Pidduck-Kent 
solution. 
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Intentionally  left  bu\nk. 
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APPENDIX  F: 

COMPUTER  PROGRAM  FOR  CALCULATION  OF 
PIDDUCK-KENT  PARAMETERS 


101 


Intentionally  left  bunk. 
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program  pidduck 
c 

c  Pidduck-Kent  constants  DEL  and  AO 
c  Note:  Accuracy  is  poor  for  charge-to-mass  ratio  <  IE-03 

c  Values  of  DEL  and  AO  are  not  calculated  for  EPS  =  0 

c  Iteration  fails  for  EPS  <  IE-05 
c  24  March  1991 

c 

c  Converted  from  BASIC  to  FORTRAN 

c  Converter:  Caledonia  Henry 

c  Date:  30  Jan  1992 

c 

dimension  a(5),  d(5) 

data  (a(i),i=1 ,5),(d(i),i=1 ,5)/-.003659837,-.C0693S81 6, 

*  .0006852334, -.00003094782, .000000529 1 406, .3667921  ,-.03838811, 

*  .003275623,  -.0001 455872, .000002498742/ 
c 

Pl=4.0*atan(1 .0) 
c 

c  Open  file  for  output 

open(9,file=’pkconst.out’) 

c 

c  Input  GAMMA  and  Charge-to-Mass  Ratio 

c 

98  write(6,,(a)')  "Input  Gamma  and  Charge-to-Mass  Ratio" 

read(5,*)  gam,  eps 
c 

c  calculate  on  and  cnn 

c 

cn  =  1.0/(gam-1.0) 
cnn  =  cn  +  1 .0 
c 

c  Calculate  initial  estimate  for  Delta  and  A0 
c  Using  Grollman  and  Baer  formula, 

c 

c  Note:  Accuracy  decreases  rapidly  for  eps  >  20 
c 

del  =  3.0 
do  100  i  =  1,5 

del  =  del+(a(i)+d(i)/gam)*(eps**i) 

100  continue 
c 

if(eps.gt.O.O)  go  to  200 
az  =0,0 
go  to  21 1 
c 

200  azl  =  (((2.0*cn)+3.0)/del)+((2.0*(cn+1.0))/eps) 
az  =  1.0/azl 
c 
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/ 


/ 


21 1  write(9,'(a)')  'GROLLMAN  AND  BAER  LEAST  SQUARES  FIT’ 
wriie(9,1000)  eps,  cn,  gam,  del,  a z 
1000  format(1x,"EPS  =  ",f8.4,*  N  =  ",{8.4,"  GAMMA  =  ",17.4,"  DEL  =  " 
*19.6,"  AO  =  ",f9,6) 
c 

if(eps.gt.O.O)  go  to  224 
azzz  =  0.0 
dell  =  3.0 
go  to  338 
c 

224  if(az.lt.0.7)  go  to  226 
a 2  =0.7 
c 

c  if  A0  >  0.65  then  use  alternate  solution 
c 

226  if(az.gt.0.65)  go  to  250 
c 

c  Newton-Raphson  iteration  for  A0  <=  0.65 
c 

228  j=0 

azz  =  az 

j  =  j+1 
c 

236  call  g0g1g2le(  g2.  azz,  cn,  eps.cnn) 
azzz  =  azz-g2 

if  (abs(g2).lt.1.E-06)  go  to  334 
c 

azz  =  azzz 

if(azz.le.0.65)  go  to  236 
az  =  azz 
c 

c  Newton-Raphson  Iteration  for  A0  >  0.65 
c 

c  First  calculate  Beta  function  -  B(0.5,CN+1.0) 
c 

250  w  =  0.5 

gu  =  sqrt(PI) 
w  =  cn  +  1 .0 
call  gamfun{gamma,w) 
gv  =  gamma 
266  w  =  cn  +  1.5 

call  gamfun(gamma.w) 
guv  =  gamma 
beta  =  gu*gv/guv 
c 

c  now  proceed  with  iteration 
c 

j=0 

azz=az 
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308  j  =  j+1 

call  g0g1g2gt  {  g2,  azz,  cn,  eps,  cnn,  beta) 
azzz  =  azz  -  g2 

if  (abs(g2)  .It.  1.E-06)  go  to  334 
azz  =  azzz 

if  (azz  .gt.  .65)  go  to  308 
az  =  azz 
go  to  228 

334  dell  =  ((1.0/azzz)-((2.0*(cn+1.0))/eps))/(2.0*cn+3.0) 
dell  =  1.0/dell 

338  write  (9/(3)')  'CALCULATION********* . 

write  (9,1000)  eps,  cn,  gam,  dell,  azzz 
c 

close(9) 

end 

c 

c 

c 

subroutine  g0g1g2le  (  g2,  azz,  cn,  eps, cnn) 
c 

c  subroutine  for  A0  <=  0.65 
c 

x=a zz 
prod  =  1 . 
fact  =  1 . 
sum  =  1. 
summ  =  0. 
k  =  0 

420  k  =  k  +  1 

prod  =  prod*(k-cn-1.) 

fact  =  fact  *  k 

arg  =  (x**k)/(2.0*k+1.0) 

argg  =  k*(x**(k-1))/(2.0*k+1 .0) 

term  =  prod*arg/fact 

termm  =  prod*argg/fact 

sum  =  sum+term 

summ  =  summ+termm 

if(abs(term).lt.1.E-10)  go  to  540 

go  to  420 

540  g0=(2.*cnn*(x/((1  ,-x)*'  v,nn))‘sum)-eps 

g1=(g0+eps)‘((1  ./x)+(cnn/(1  ,-x)))+2.*crm*(x/((1 .0-x)‘*cnn))*summ 
g2  =  go/gi 
return 
end 
c 
c 
c 

subroutine  gOg  1  g2gt  (  g2,  azz,  cn,  eps,  cnn,  beta) 
c 
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c  subroutine  for  AO  >  0.65 
c 

x  =  1 .0-azz 
prod  =  1 . 
fact  =  1 . 
sum  =  1. 
summ  =  0. 
k  =  0 

670  k  =  k  +  1 

prod  =  prod*(k+0.5-1) 
fact  =  fact'k 

arg  =  (x**k)*(cnn/(k+cnn)) 
argg  =  k*(x**(k-1))*{cnn/(k+cnn)) 
term  =  prod*arg/fact 
termm  =  prod*argg/fact 
sum  =  sum+term 
summ  =  summ-termm 
if(abs(term).lt.1.E-10)  go  to  800 
go  to  670 

800  ssum  =  (0.5/azz"0.5)*(beta-{{x“cnn)/cnn)*sum) 

ssumm  =  ((-0.5/azz)*ssum)+{(0.5/azz**0.5)*(x**cn)*(sum-{x/cnn)‘ 

*  *summ)) 

gO  =  (2.0*cnn*(azz/(x**cnn))*ssum)-eps 
gl  =  (g0+eps)*({1.0/azz)+(cnn/x))+2.0*cnn*{azz/(x**cnn))*ssumm 
g2  =  gO/gl 
return 
end 
c 
c 
c 

subroutine  gamfun(gamma,w) 
c 

c  gamma  function 
c 

dimension  b(8) 
c 

data  (b{i),i=1 ,8)/-.5771 9 1 652,  .988205891,  -.897056937, 
*.918206857,-756704078,  .482199394,  -.193527818,  .035868343/ 
c 

iw  =  int{w) 
fw  =  w-iw 

if(iw.gt.O)  go  to  920 
gamma  =  1 .0/fw 
go  to  970 
920  gamma  =1 

do  960  j  =  1,iw-1 
gamma  =  gamma'(w-j) 

960  continue 
970  gfw  =  1 
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do  1000  k=  1,8 
gfw  =  gfw  +  (b{k)*(fw**k)) 
1000  continue 

gamma  =  gamma'gfw 

return 

end 


107 


INTENTIONALLY  LEFT  BLANK. 


LIST  OF  SYMBOLS 


ao 

B 

Bx(a,b) 

C 

C(x0) 

<V 

Ej(t) 

C 

^0 

2F1(a,b;c;z) 

f(x0) 

fb  s  f(v0) 

K,  K(x0) 

KE 

M 

n 

P 

W) 

P(x0.  t) 

q(p) 

R 

t 

T 

W) 

T(x0.t) 

u(x0.t) 

UP 

W(x0,t) 

x(x0,t) 

xo 

y(t)  -  x(y0,t) 
Vo 

z(x0,t) 

y 


Bore  Area 

Parameter  (Eq.  B-28) 

Pidduck-Kent  Characteristic  Parameter  (Eq.  B-33) 

Separation  Constant  (Eq.  B-23) 

Incomplete  Beta  Function 

Propellant  Charge  Mass 

Gas  Mass  Between  Breech  and  x0  (Eq,  B-113) 

Specific  Heat  at  Constant  Volume 
Gas  Internal  Energy 
Initial  Gas  Internal  Energy 
Gauss  Hypergeometric  Function 

Function  Describing  Spatial  Dependence  of  z  (x0,t)  (Eq.  B-19) 

Value  of  f(x0)  at  Projectile  Base 

Adiabatic  Constant 

Kinetic  Energy 

Projectile  Mass 

Polytropic  Index,  (7- I)"1 

Gas  Pressure 

Space  Mean  Gas  Pressure 

Pressure  of  Gas  Element  Initially  Located  at  x0 

integration  Variable  (Eq.  B-41) 

Gas  Constant 

Integral  Arising  in  Pidduck-Kent  Derivation;  Specific  Gauss  Hypergeometric 
Function  (Eq.  B-40  and  B-41) 

Time 

Gas  Temperature 

Space  Mean  Gas  Temperature 

Temperature  of  Gas  Element  Initially  Located  at  x0 

Velocity  of  Gas  Element  Initially  Located  at  x0 

Projectile  Velocity 

Function  of  Density  of  Gas  Element  Initially  Located  at  x0  (Eq.  B-12) 

Position  of  Gas  Element  Initially  Located  at  x0 

Initial  Location  of  Gas  Element 

Projectile  Position 

Initial  Projectile  Position 

Pidduck-Kent  Variable,  Reduced  Gas  Element  Position  (Eq.  B-7) 

Ratio  of  Specific  Heats 
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5  Pidduck-Kent  Characteristic  Parameter  (Eq.  B-83) 

e  Ratio  of  Propellant  Charge  Mass  to  Projectile  Mass 

r|  Gas  Ccvolume 

p(x0)  Integration  Variable  (Eq.  B-68) 

p  Gas  Density 

p(x0,t)  Density  of  Gas  Element  Initially  Located  at  x0 

p(f)  Space  Mean  Gas  Density 
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3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of 
ideas,  etc.)  _ _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 
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U.S.  Army  Research  Laboratory  paBmaa 

ATTN:  AMSRL-OP-CI-B  (Tech  Lib) 

Aberdeen  Proving  Ground,  MD  21005-5066 
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